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In this work the determination of low-energy bound states in Quantum Chromodynamics is recast so that it 

is linked to a weak-coupling problem. This allows one to approach the solution with the same techniques which 

solve Quantum Electrodynamics: namely, a combination of weak-coupling diagrams and many-body quantum 

mechanics. The key to eliminating necessarily nonperturbative effects is the use of a bare Hamiltonian in which 

quarks and gluons have nonzero constituent masses rather than the zero masses of the current picture. The use 

^^ ' of constituent masses cuts off the growth of the running coupling constant and makes it possible that the running 

Q^ ^ coupling never leaves the perturbative domain. For stabilization purposes an artificial potential is added to the 

^\ . Hamiltonian, but with a coefficient that vanishes at the physical value of the coupling constant. The weak-coupling 

approach potentially reconciles the simplicity of the Constituent Quark Model with the complexities of Quantum 

(— I , Chromodynamics. The penalty for achieving this perturbative picture is the necessity of formulating the dynamics 

j^ ' of QCD in light-front coordinates and of dealing with the complexities of renormalization which such a formulation 

entails. We describe the renormalization process first using a qualitative phase space cell analysis, and we then 

QQ , set up a precise similarity renormalization scheme with cutoffs on constituent momenta and exhibit calculations to 

^SJ ' second order. We outline further computations that remain to be carried out. There is an initial nonperturbative 

but nonrelativistic calculation of the hadronic masses that determines the artificial potential, with binding energies 

required to be fourth order in the coupling as in QED. Next there is a calculation of the leading radiative corrections 

^ ' to these masses, which requires our renormalization program. Then the real struggle of finding the right extensions 

rn : to perturbation theory to study the strong-coupling behavior of bound states can begin. 
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I. INTRODUCTION 



O ; The only truly successful approach to bound states in field theory has been Quantum Electrodynamics (QED), 

^ ' with its combination of nonrelativistic quantum mechanics to handle bound states and perturbation theory to handle 
,.^ relativistic effects. Lattice Gauge Theory is maturing but has yet to rival QED's comprehensive success. There 

are four barriers which prohibit an approach to Quantum Chromodynamics (QCD) that is analogous to QED. The 
barriers are: (1) the unlimited growth of the running coupling constant g in the infrared region, which invalidates 
perturbation theory; (2) confinement, which requires potentials that diverge at long distances as opposed to the 
Coulombic potentials of perturbation theory; (3) spontaneous chiral symmetry breaking, which does not occur in 
perturbation theory; and (4) the nonperturbative structure of the QCD vacuum. Contrasting the gloomy picture of 
the strong interaction in QCD, however, is that of the Constituent Quark Model (CQM), where only the minimum 
number of constituents required by the symmetries are used to build each hadron and where Zweig's rule leaves little 
role for production of extra constituents. Instead, rearrangement of pre-existing constituents dominates the physics. 
Yet the CQM has never been reconciled with QCD — not even qualitatively. Ever since the work of Feynman, though, 
it has been clear that the best hope of reconciliation is offered by infinite momentum frame (IMF) dynamics. 

In this paper, a framework closely related to the IMF will be employed, namely, light-front quantization. The 
purpose of this paper is to provide arguments that all the barriers to a perturbative starting point for solving QCD 
can be overcome, at least in principle, when a light-front framework is used. We present the basic formulation for 
such an approach. Coupled to the choice of coordinates are the introduction of nonzero masses for both quarks and 
gluons, the use of cutoffs on constituent momenta which eliminate vacuum degrees of freedom, and the addition of 
an artificial stabilizing and confining potential which vanishes at the relativistic value of the renormalized coupling 
but nowhere else. One result of these unconventional modifications is a theory with a trivial vacuum. However, the 
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theory has extra terms in the Hamiltonian which are induced by the ehmination of vacuum degrees of freedom, which 
account for spontaneous chiral symmetry breaking, yet can be treated perturbatively. It is a theory where a second- 
order treatment of renormahzation effects should closely resemble the phenomenological CQM, while a fourth-order 
treatment — if all goes well — should begin to replace phenomenology by true results of QCD. We have not carried 
out these computations but will outline them to show that they are indeed perturbative ones. 

Our basic aim, then, is to tailor an approach to QCD which is based on the phenomenologically successful approach 
of the CQM. Central to our formulation is to start with nonzero masses for both quarks and gluons and to then 
consider the case of an arbitrary coupling constant g which is small even at the quark-gluon mass scale — a running 
coupling which is then small everywhere because below this scale it cannot run anymore. This means sacrificing 
manifest gauge invariance and Lorentz covariance, with these symmetries only being implicitly restored (if at all) 
when the renormalized coupling is increased to its relativistic value, which we call g^. The value gs is a fixed number 
because it is measured at the hadron mass scale, and by asymptotic freedom arguments such a coupling has a fixed 
value. For smaller values of g our theory lacks full covariance and is not expected to have the predictive power of 
QCD, but it allows phenomenology to guide renormahzation and is defined to maximize the ease of perturbative 
computations and extrapolation to gs- 

Some key new ideas in this paper have previously been reported only in unpublished notes m^- Some have been 
presented in a condensed form in Ref. q]. These ideas underly two main parts of this paper. The first part is a 
qualitative power counting analysis of divergences in light-front QCD (LFQCD), which provides a vital basis for the 
new cutoff scheme and renormahzation framework we develop in the second part. More broadly, this paper draws on 
a range of previous research efforts by the authors and colleagues: the QCD calculations presented in this paper are 
largely based on the standard gauge-fixed LFQCD Hamiltonian |Q, and the specific diagrammatic rules used here are 
defined with examples in Refs. [pj^f^ following earlier work [lO]. This paper follows recent work on QCD |P-[l^,0,0| 



and the similarity renormahzation scheme [114 15 1; and it is also the outgrowth of a line of research that covers the 
Tamm-Dancoff approximation [|l6|, light-front QED [ll^, and light-front renormahzation ||l^. The relation of this 
paper to the entire field of light-front field theory will be outlined in section II. 

The plan of the paper is as follows. As our formulation is quite different not only from previous approaches 
to QCD but also to previous studies of light-front dynamics, we start in Section II with a general overview and 
motivation of our approach. We discuss the apparent contradictions between QCD and its precursor, the CQM. Then 
we outline why a light-front approach can bridge the gap between QCD and the CQM and discuss in turn how each 
of the equal-time barriers to a weak-coupling treatment of QCD may be overcome. In Sections III-V, the first major 
part of the paper, we motivate and discuss the light-front power-counting analysis, describe the canonical LFQCD 
Hamiltonian, and use a qualitative phase space cell analysis to classify the possible counterterms for ultraviolet and 
infrared divergences in LFQCD needed for renormahzation. Light-front power counting can determine the operator 
structure of the LFQCD Hamiltonian 0, and we rely on the derivation from the QCD Lagrangian to determine 
details such as color factors. This canonical LFQCD Hamiltonian will be treated as a Hamiltonian in the many-body 
space of finite numbers of quark and gluon constituents, a Hamiltonian which with regularization gives finite answers 
in this many-body space. Three terms in the canonical Hamiltonian itself require infrared counterterms because of 
linear infrared divergences, and we proceed to construct them. We explain the absence of finite parts to many of these 
counterterms because kinematical longitudinal boost invariance is a scale invariance. Because of this kinematical 
constraint, only counterterms for logarithmic infrared divergences are allowed to contain finite parts, which as part of 
the renormahzation process must be tuned to reproduce physical results. Counterterms to canonical terms are prime 
candidates for bringing in phenomena associated with true confinement and with the spontaneously broken vacuum 
of normal rest frames. We distinguish true confinement in the exact theory from the artificial confinement which we 
introduce by hand in the weak-coupling starting point. 

This leads to the second major part of the paper, where in Sections VI-X we set up an explicit quantitative scheme 
for renormahzation of the LFQCD Hamiltonian. First, a momentum space cutoff procedure is introduced which 
regulates the Hamiltonian itself rather than the individual terms of perturbation theory Q. This cutoff scheme 
depends on the use of massive constituents and is chosen to ensure that each state has only a finite number of 
constituents. Both gauge invariance and Lorentz covariance are violated by these cutoffs, and a range of counterterms 
is needed to enable a finite limit as the cutoffs are removed. Restoration of the violated symmetries can only be 
established by examining the solution to the theory at the relativistic limit gs (or good approximations to this 
solution). This is the price we must pay for achieving the new framework. We proceed to the construction of the 
effective Hamiltonian. We begin by discussing a novel perturbation theory formalism due to Glazek and Wilson 
||l4| , [l5| designed to transform a cutoff Hamiltonian to band-diagonal form while avoiding small energy denominators 
that plague other approaches |l9|. This is called the "similarity renormahzation scheme." The end product of the 
scheme is a band-diagonal effective Hamiltonian in which dependence on the original cutoffs has been removed to any 
desired order of perturbation theory. As an example we determine in detail the gluon mass counterterms necessary 
to remove second-order divergences generated by two-gluon intermediate states. We then proceed to construct order 



g^ light-front infrared counterterms, and we identify finite counterterms that may be necessary to compensate for 
the removal of zero longitudinal momentum modes from the cutoff theory. The next step is to analyze the role of an 
artificial potential in the model. We use it to maintain the qualitative structure of bound states for small g, except 
for an overall scale factor g'^. We argue our theory is closely analogous to a CQM if the band-diagonal Hamiltonian 
is computed only to order g^ . Our hope is that higher-order computations for the effective Hamiltonian will be 
invaluable preparation for study of the limit g ^f ga- 
ll. MOTIVATION AND OVERVIEW 

A key ingredient in the constituent picture of hadrons, starting with the parton model and then moving to the quark 
parton model with constituent quarks and gluons, is the infinite momentum frame. The field theoretic realization of 
the intuitive ideas from the IMF is provided by the light-front dynamics of field theories. Recently there has been 
renewed interest in LFQCD because of its potential advantages over the normal equal-time formulation, especially 
due to the triviality of the light-front vacuum in the cutoff theory. However, a basic puzzle has remained — namely, 
how do confinement, chiral symmetry breaking, and other nonperturbative as pec ts of QCD emerge from LFQCD? 

Since Dirac's formulation of Hamiltonian systems in light-front coordinates pS and the development of the infinite 



momentum frame limit of equal-time field theory | 
by the recognition of its importance for current a 



21| , there has been intermittent progress in this area, initially driven 
gebra |2^ and the parton model |2^ , but later slowed by the many 
renormalization problems inherent to Hamiltonian field theory. After early progress towards understanding light-front 
field theory ||2j], theorists began to consider light-front QCD |0|. No clearly successful approach to nonperturbative 
QCD emerged, so some theorists tried to make progress on technically simpler problems such as 1-1-1 dimensional 
field theory |2^] and bound states in QED and the Yukawa model in 3-1-1 dimensions |26|. More recently there has 
also been much work on the light-front zero mode problem ||27|] , primarily by theorists who advocate a different 
approach to the vacuum problem than that developed in this paper — the relationship between most of this work and 
ours is at present obscure. Unfortunately, the barriers that are discussed in this paper have hindered theorists from 
attacking nonperturbative QCD in 3-1-1 dimensions; and so recent work on light-front field theory has not focused on 
the light-front problem that is of greatest practical interest. 

In this paper we establish a new framework for studying QCD in light-front coordinates by building from an 
unconventional choice of the bare cutoff Hamiltonian. The basic question we try to answer is the following: can 
one set up QCD to be renormalized and solved by the same techniques that solve QED; namely, a combination of 
weak-coupling perturbation theory and many-body quantum mechanics? The changes from the standard approach in 
QED we introduce are: a) the use of light-front dynamics; b) the use of nonzero masses for both quarks and gluons; 
c) the need to handle relativistic effects which give rise to, for example, asymptotic freedom in QCD, which in turn 
leads to a fairly strong renormalized coupling constant and hence relativistic binding energies; d) the presence of 
artificial stabilizing and confining potentials which vanish at relativistic values of the coupling constant but nowhere 
else; and e) the concerns about light- front longitudinal infrared divergences, which cancel perturbatively in QED 
because of gauge invariance. Our theory is not gauge-invariant order by order because we use a non-zero gluon mass. 
Moreover, we propose that these non-cancelling divergences are necessarily the sources of true confining potentials 
and chiral symmetry breaking in QCD. With the cutoff Hamiltonian we have a trivial vacuum. Furthermore, the 
free part of the cutoff Hamiltonian exhibits exact chiral symmetry despite the existence of a quark mass. (In the 
canonical light-front Hamiltonian, chiral symmetry is explicitly broken only by the helicity flip part of the quark- 
gluon interaction. The free Hamiltonian itself only breaks chiral symmetry if zero mode quarks are included.) We 
expect light-front infrared divergences to be sources of confinement and chiral symmetry breaking because these are 
both vacuum effects in QCD, and we show that vacuum effects can enter the light-front theory through light-front 
longitudinal infrared effects. Because of the unconventional scaling properties of the light-front Hamiltonian, these 
effects include renormalization counterterms with whole functions to be determined by the renormalization process. 

The basic motivation for our approach is physical rather than mathematical. Physically, one's unperturbed or 
starting Hamiltonian is supposed to model the physics one is after, at least roughly. A Hamiltonian with nonzero 
(constituent) quark and gluon masses and confining potentials is closer to the physics of strong interactions than 
a Hamiltonian with zero mass constituents and no confining potentials. Then, in the spirit of QED, we analyze 
renormalization effects with the confining potential itself treated perturbatively, but only to generate an effective 
few-body Hamiltonian which can be solved nonperturbatively. 



A. Constituent quarks versus QCD 

Prior to the establishment of QCD as the underlying theory of strong interactions, there arose the Constituent Quark 
Model Ig^ and Feynman's Parton Model ^^. The CQM provides an intuitive understanding of many low-energy 
observables. The Parton Model provides an intuitive understanding of many high-energy phenomena. Then QCD 
came along, and along with it came the commonly accepted notion that the vacuum of QCD is a very complicated 
medium. Unfortunately, this is the source of a contradiction between the constituent quark picture and QCD; and 
twenty years of study of QCD have done little to weaken that contradiction. 

The complicated vacuum of QCD plays a crucial role in invalidating any perturbative picture of isolated quarks and 
gluons. One important function of the vacuum is to produce confining interactions among quarks and gluons at large 
distances, thus overturning the non-confining gluon exchange potential of perturbation theory. Another important 
function of the vacuum is to provide the spontaneous breaking of chiral symmetry. The equal-time QCD vacuum is an 
infinite sea of quarks and gluons, and baryons and mesons arise as excitations on this sea. Unfortunately, individual 
quarks and gluons are lost in this infinite sea. 

According to the CQM, however, a meson is a simple quark-antiquark bound state and a baryon is a bound state 
of three quarks; and Zweig's rule suppresses particle production in favor of rearrangement of constituents. But how 
can the hadrons be simply a quark-antiquark or a three quark bound state if they are excitations over a complicated 
vacuum state? One may hope that the effects of quarks and gluons in the vacuum may be treated via weak-coupling 
perturbation theory. Unfortunately, the coupling of quarks and gluons in the vacuum grows in strength as the average 
momentum of these constituents decreases as a result of asymptotic freedom. Thus one expects the density of low- 
momentum constituents in the vacuum to be quite large, thereby invalidating any perturbative treatment of vacuum 
effects. 

The CQM, as we conceive it, requires that both quarks and gluons have sizable masses. For gluons, this violates 
the gauge invariance of QCD. For quarks, this violates the rule that chiral symmetry is not explicitly broken only for 
massless quarks. 

On the other hand, many high-energy phenomena are most naturally described in the language of the Parton Model. 
The constituent picture and the probabilistic interpretation of distribution functions are essential for the validity of 
the Parton Model. But it is not at all easy to reconcile the probabilistic picture with the notion of a nontrivial vacuum 
in the equal-time framework. Thus both the Constituent Quark Model and the Parton Model are put in peril by 
QCD with a complicated vacuum structure. 

At first glance, the blame for the contradiction rests with the naivete of the constituent picture. Equal-time QCD 
seems to clearly indicate that any few-body, perturbative approach to hadron bound states is unfounded. However, 
in light-front dynamics we argue that the CQM and QCD can be reconciled so that the apparent contradictions 
disappear. 

B. Why the Ught front and massive constituents? 

A constituent picture of hadrons is certainly very natural in a nonrelativistic context. However, particle creation 
and destruction need to be a part of any realistic picture of relativistic bound states. Is it possible to build a relativistic 
constituent picture of hadrons based on nonrelativistic, few-body intuitions? 

As pointed out long ago, the most serious obstacles to this goal are overcome by changing to light-front coordinates 
(the "front-form" in Dirac's original work POJ), or moving to the IMF |23,p4,El|. In light-front quantization one 
quantizes on a surface at fixed light-front time, x'^ = t + z (see Fig. 1), and evolves the system using a light-front 
Hamiltonian P^, which is the momentum conjugate to a;+. A longitudinal spatial coordinate, x^ = t~ z, arises, with 
its conjugate longitudinal momentum being P^. 

The possibility of building a constituent picture in light-front field theory rests on a simple observation. All physical 
trajectories lie in or on the forward light cone. This means that all physical trajectories lie in the first quadrant of 
the light-front coordinate system, so that all longitudinal momenta satisfy the constraint, 

fc+ > 0. (2.1) 

By implementing any cutoff that removes degrees of freedom with identically zero longitudinal momentum, one forces 
the vacuum to be trivial because it can carry no longitudinal momentum. 
For a free massive particle on shell (A;^ = m^), the light-front energy is 

,- ^ fi±=^. (2.2, 



where k± is the transverse momentum. This means that the zero-momentum states we must remove to create a trivial 
vacuum in theories with positive m?, have infinite energy unless k± = and m = 0. This makes it sensible to replace 
the zero-momentum modes with effective interactions, since this is exactly the strategy used when renormalizing 
divergences from high energy degrees of freedom in equal-time field theory. However, such a starting point may be 
far removed from the canonical field theory. 

When field theories simpler than QCD are analyzed in light-front coordinates, it becomes apparent that the as- 
sumption of a trivial vacuum can be misleading. If m^ < because of spontaneous or dynamical symmetry breaking, 
constituents may have exactly zero longitudinal momentum (known as "zero modes" ) and still have finite energy. 
In this case the physical vacuum is free to contain an arbitrarily large number of zero longitudinal momentum con- 
stituents. The importance of zero modes is most simply illustrated in models with spontaneous symmetry breaking 
such as the sigma model. In the case of the sigma model, however, using power-counting arguments and demanding 
current conservation (see Appendix A) one can easily determine the counterterms needed to preserve the physics — 
at least on the canonical or tree level — for a theory in which zero modes are removed. We compare and contrast the 
sigma model with and without the zero modes removed. Both are reasonable theories, but the phenomena of the two 
theories must be described with different languages. The discussion of Appendix A indicates that an alternative way 
to realize the dynamics of spontaneous symmetry breaking on the light-front is to force the vacuum to be trivial and 
to include counterterms in the Hamiltonian which are based on power counting and explicitly break chiral symmetry. 
In Appendix A it is shown that current conservation can be used to fix these counterterms. It has also been shown in 
Rcf. [ p2| that in some cases one can use the renormalization group and coupling coherence, which is discussed below, 
to fix such counterterms. 

The analysis of QCD is far more complicated than that of the sigma model. In QCD, the signal of spontaneous 
symmetry breaking is still a vanishing pion mass, but the pion is now a composite state. For the weak-coupling 
starting point, we require all hadrons including the pion to have masses that are close to the sum of their constituent 
quark masses. One might imagine other limiting procedures, but we require our starting point to be perturbative 
and the pion cannot be massless initially with this criterion. As the coupling increases, the pion mass must decrease 
toward zero if spontaneous chiral symmetry breaking is to be recovered ss g —^ gs- The pion mass should be a 
continuous function of the coupling, so it cannot reach zero until the coupling reaches its physical value gs] and it does 
so only after the inclusion of renormalization effects. In contrast, the sigma model illustrates spontaneous symmetry 
breaking even for arbitrarily small coupling. Thus the complete determination of the terms necessary to counter the 
elimination of zero modes in QCD will not be simple. 

But the major problem in LFQCD is not the question of zero modes. To even address the role of zero modes we need 
a reliable, practical calculational framework. LFQCD has severe infrared divergences arising from small longitudinal 
momenta when we eliminate the zero modes [H,E2lH-nOJ . These infrared divergences are separate from and in addition 
to the infrared problems of equal-time QCD. In equal-time QCD, infrared problems arise due to both zero quark and 
gluon masses and to the growth of the running coupling constant in the infrared domain. In LFQCD the same infrared 
problems also exist, but they are divergences associated with small transverse momenta, which are the only momenta 
that combine directly with masses. The longitudinal infrared divergences are special to LFQCD and for this reason it 
is normally presumed that they will cancel out if treated properly. An essential part of this supposed cancellation is 
the maintenance of gauge invariance. To preserve manifest gauge invariance in QCD in perturbation theory one needs 
massless gluons and carefully chosen regularization schemes. With massless gluons, however, the running coupling 
constant increases without bound when the energy scale is of the order of hadronic bound state energies. Thus one 
needs all orders of perturbation theory to compute observables in the hadronic bound state range. But all orders of 
perturbation theory involve arbitrary numbers of quarks and gluons as intermediate states, thus contradicting the 
notion that hadrons are mostly a quark-antiquark or three-quark bound state. It is this problem that has caused us 
to abandon manifest gauge invariance in favor of a weak-coupling picture in which gluons have mass. 

Suppose we consider a light-front Hamiltonian whose free part corresponds to that of massive quarks and gluons. 
What is the justification for taking massive gluons in the free part of the Hamiltonian? As is always the case, the 
division of the Hamiltonian into a free part and an interaction part is arbitrary; however, it is also true that the 
convergence of a perturbative expansion depends crucially on how this choice is made. We are free to take advantage 
of this arbitrariness; we choose the free gluon part of the bare cutoff Hamiltonian to be that of massive gluons. (In 
QED, in contrast, where electrons and photons appear as free particles in asymptotic scattering states, it would be 
more difficult to exercise this freedom.) Furthermore, as we show in Section XIII, the cutoffs we employ inevitably 
introduce quadratic mass renormalization for quarks and gluons; thus the bare quark and gluon masses are tunable 
parameters. Of course, the crucial question is what the renormalized masses used for the bound state computations 
will be: these might be of the order that phenomenology assigns to constituent quarks and gluons. (In QED, one 
would tune the bare masses to reproduce physical masses; in QCD, we must tune to fit bound state properties.) 

The fact that we are giving the gluon a mass should not create any contradiction with asymptotic freedom when g 
achieves its relativistic value gg- The reason is that gs is a running coupling constant, gA, at a scale where A is of the 



same magnitude as hadronic masses, A ^^ Aqcd- 9 a is small at extremely large momentum scales, and the running 
scale is A « Aqcds'^ '^ for small running g^, where c is a positive constant. But changes to amplitudes due to masses 
can be treated perturbatively at such scales and behave as powers of Aqcd/-^ ~ e~'^l^'^, which vanishes to all orders 
in a perturbative expansion in powers of ^a- 

Once we assume the free part to consist of massive gluons, what are the consequences? A gluon mass automatically 
prevents unbounded growth of the running coupling constant below the gluon mass scale and provides kinematic 
barriers to unlimited gluon emission. It eliminates any equal-time type infrared problems. With nonzero quark and 
gluon masses it is also possible to develop a cutoff procedure for the Hamiltonian such that if the cutoffs are imposed 
in a specific frame, a large number of states (the upper limit of whose invariant masses are guaranteed to be above 
a large cutoff) are still available for study even in boosted frames. Another consequence of nonzero gluon mass is 
that the long-range gluon exchange potential between a pair of quarks is too small in transverse directions, falling off 
exponentially. Hence an artificial potential must be added to provide a full strength potential to yield realistic bound 
states for small g. 

C. Light- front infrared divergences 

Our basic objective is to establish a weak-coupling framework for studying QCD bound states, so that one can 
smoothly approach the strong-coupling limit and use bound state phenomenology to guide renormalization. We have 
argued that the choice of light-front dynamics, massive quarks and gluons, and a particular cutoff scheme eliminates 
the traditional barriers to a perturbative treatment of QCD. The first step, then, is the construction of a bare 
Hamiltonian which incorporates confining potentials, massive quarks and gluons, and a trivial vacuum. The cutoffs 
will violate Lorentz and gauge symmetries, forcing the bare Hamiltonian to contain a larger than normal suite of 
counterterms to enable a finite limit as the cutoffs are removed. 

Now in the equal-time theory, the QCD vacuum is thought to be a complicated medium which presumably provides 
both confinement and the spontaneous breaking of chiral symmetry. But in the light-front theory with a suitably 
cutoff Hamiltonian, we have asserted that the vacuum is trivial. So we have to find other sources for confinement and 
spontaneous chiral symmetry breaking in the cutoff theory. A natural place to look for them is in the divergences 
associated with light-front infrared singularities. 

Explicitly, suppose we regulate infrared divergences (where fc,^ — + 0) by cutting off the longitudinal momentum of 
each constituent i so that fc^ > e, with e a small but finite positive constant. (We will also need to regulate ultraviolet 
divergences by cutting off large transverse momenta, for example via k\ < A^ , but this is not important for the present 
discussion.) Because the total longitudinal momentum P+ of a state is just the sum of the longitudinal momenta 
of its constituents, P+ — X^i^i^j if follows immediately that the vacuum of the cutoff theory (for which P+ = 0) 
contains no constituents; that is, the vacuum is trivial. Moreover, any state with finite P^ can contain at most 
P"'"/e constituents. Therefore, effects which in an equal-time formulation are due to infinite numbers of constituents 
— in particular, confinement and chiral symmetry breaking — must have other sources in the cutoff theory on the 
light-front. The obvious candidates are the counterterms which must be introduced in the effective Hamiltonian in 
order to eliminate the dependence of observables on the infrared cutoff e. Of special interest are counterterms that 
reflect consequences of zero modes (namely, modes with fc+ — 0) in the full theory (see discussions in Sections IV. B 
and XI). 

Now in the canonical Hamiltonian, one particular term of interest is the instantaneous interaction in the longitudinal 
direction between color charge densities, which provides a potential which is linear in the longitudinal separation 
between two constituents that have the same transverse positions. In the absence of a gluon mass term, this interaction 
is precisely cancelled by the emission and absorption of longitudinal gluons. The presence of a gluon mass term means 
that this cancellation becomes incomplete. But what is still lacking is the source of transverse confining interactions. 

According to power-counting arguments, the counterterms for longitudinal light-front infrared divergences may 
contain functions of transverse momenta B; and there exists the possibility that the a priori unknown functions 
in the finite parts of these counterterms will include confining interactions in the transverse direction. The g — > gs 
limit may then be smooth if such confining functions are actually required to restore full covariance to the theory. In 
the absence of a gluon mass term, the light-front singularities are supposed to cancel among each other in physical 
amplitudes. This has been verified to order g^ explicitly in perturbative amplitudes for quarks and gluons rather than 
physical states, to order g'^ for the quark-gluon vertex |l2[[lC| ], and to order g'^ in the gluon four-point function [p3[ . 
A gluon mass term in the free part of the Hamiltonian, however, spoils this cancellation. What are the consequences? 

The light-front infrared singularities give rise to both linear and logarithmic divergences. The linear divergences, 
however, contain the inverse of the longitudinal cutoff i, which violates longitudinal boost invariance; and hence the 
infinite parts of the counterterms for these divergences also violate longitudinal boost invariance. Thus finite parts 



for these counterterms are prohibited by longitudinal boost invariance, which is a kinematical symmetry. So we have 
to analyze logarithmic infrared divergences in order to get candidates for transverse confinement. 

There is a specific problem that the complexity of the counterterms creates. In covariant perturbation theory 
every coupling or mass introduced by renormalization becomes an independent parameter. One might worry that the 
appearance of whole functions as counterterms could destroy the predictive power of the theory because functions 
include an infinite number of parameters and may seem to destroy the renormalizability of the theory. We discuss 
the resolution of this problem at the end of Section V. 

D. Zero modes and chiral symmetry breaking 

How does chiral symmetry breaking arise in the light-front theory with the zero modes removed? The removal of 
the zero modes has two important consequences for chiral symmetry in LFQCD — namely, the vacuum is trivial, and 
chiral symmetry is exact for free quarks of any mass. The consequence of a trivial vacuum is that, as with the sigma 
model, the mechanism for effects associated with spontaneous chiral symmetry breaking in the equal-time theory 
will be far different in LFQCD. The consequence of chiral symmetry being exact for massive constituents means 
that the mechanism for effects associated with explicit chiral symmetry breaking in equal-time will also be quite 
different in LFQCD. We discuss these points further in this subsection but leave most details of the light-front chiral 
transformation to Appendix B. We just note here for the discussion which follows that the fcrmion field naturally 
separates into two-component fields ip — "0+ -I- ?/;_ on the light-front, where "0+ is dynamical and ip- is constrained. 
The light-front chiral transformation applies freely only to the two-component field tp^ |34|, because the constraint 
equations are inconsistent with the chiral transformation rules when explicit breaking is present. 

In light-front dynamics chiral symmetry is exact for free quarks of any mass once zero modes have been removed. 
This is because chiral charge conservation is simply the conservation of light-front helicity, which is a fundamental 
property of free quarks in light front dynamics in the absence of zero modes (see Appendix B). Despite the conservation 
of chiral charge, the local chiral current is not conserved. The divergence of the chiral current remains 2impil)"f^'ip. 
How does one reconcile the conservation of light-front chiral charge with the non-conservation of the axial current? 
From the chiral current divergence equation, the light-front time derivative of the chiral charge is proportional to 
J dx~(Px±'ijj^5Tp. When the fields in this integral are expanded in terms of momentum eigenstates, the diagonal terms 
b*b and d*d — where b and d are the quark and antiquark annihilation operators — vanish, namely the matrix elements 
multiplying them vanish. Moreover, the off-diagonal terms b^d^ , and bd vanish if the zero modes are absent. Thus it 
is the absence of zero modes which makes it possible for the light-front chiral charge to be conserved irrespective of 
the non-conservation of the axial current for free massive fermions. The light-front time derivative of the chiral charge 
can avoid vanishing only if "zero mode quarks and antiquarks" (quarks and antiquarks with exactly zero longitudinal 
momentum) are permitted to exist. But the cutoffs we use prevent this possibility, and hence chiral symmetry is 
exact for all free quarks in the cutoff theory. 

In normal reference frames, the absence of a mass term implies conservation of the axial vector current and hence 
a conserved axial charge. Given a conserved axial current, there are two possibilities: a) the axial charge annihilates 
the vacuum, in which case as a consequence of Coleman's theorem ( "invariance of the vacuum is the invariance of the 
world" ) the symmetry is reflected in the spectrum of the Hamiltonian — that is, we expect degenerate parity doublets 
in the spectrum; or b) the axial charge does not annihilate the vacuum, in which case one talks about the spontaneous 
breaking of chiral symmetry, as a consequence of which massless Nambu-Goldstone bosons should exist. In the real 
world, the pion — the Goldstone boson — is very light, and the second possibility is thought to be realized. The 
nonzero mass of the pion is thought to arise from the explicit symmetry breaking terms in the Lagrangian, namely, 
the small current quark masses. 

On the light-front, in the absence of interactions, massive quarks in the cutoff theory do not violate the chiral 
symmetry of the light-front Hamiltonian thanks to the cutoffs. Now consider interactions in the gauge theory. There 
remains one explicit chiral breaking term in the canonical Hamiltonian (see Section IV), 



gmp / dx d^x^il)\(7±_ ■ iA^—i}}+ - — (A_lV+) j . 



(2.3) 



This term, which involves gluon emission and absorption, is linear in the quark mass and couples the transverse gluon 
field A± to the Dirac matrix 7j_ (actually, the Pauli matrix a± in the two-component notation above), which causes 
a helicity flip. In the free quark Hamiltonian, 7j^ does not appear, and the quark mass appears only squared. Now 
the bare cutoff Hamiltonian has canonical terms and counterterms. The chiral charge still annihilates the vacuum 
state, which is just a kinematical property of the cutoff theory, despite the fact that the chiral charge no longer 
commutes with even the canonical cutoff Hamiltonian. Thus the vacuum annihilation property of the chiral charge in 



the cutoff theory has nothing to do with the symmetry of the Hamiltonian and does not appear to have any dynamical 
consequences. The concept of spontaneous symmetry breaking seems to have lost its relevance in the cutoff theory. 

But the pion should emerge as an almost massless particle. How does this become possible on the light-front without 
zero modes? We may take a hint from the sigma model with the zero modes removed, as discussed in Appendix A. In 
that model terms which explicitly violate the symmetry must be added to the Hamiltonian to yield a conserved chiral 
current, and at the same time the pion must be held massless. In QCD, the elimination of these modes directly results 
in effective interactions. The effective interactions (the counterterms in the cutoff theory) can explicitly violate chiral 
symmetry yet still persist in the limit of zero quark mass. Since spontaneous breaking of chiral symmetry in normal 
frames is a vacuum effect, we look for interaction terms that are sensitive to zero mode quarks and antiquarks. There 
is a term in the canonical Hamiltonian of the form 



g' 



/ dx (fx^iPla^ ■ A^ — (o-i • Ai V'+) , (2.4) 



which contains an instantaneous fermion. Some of the interactions which correspond to this term arc shown dia- 
grammatically in Fig. 2. These interactions by themselves do not violate chiral symmetry. However, since they are 
sensitive to fermion zero modes, the counterterms for these interactions need not be restricted by chiral symmetry 
considerations. The only requirement we can impose is that they obey light-front power counting criteria. 

We show later that there are explicit chiral symmetry breaking terms which satisfy the power counting restriction. 
We add such terms to our Hamiltonian. As a result of renormalization, in addition to the emergence of non-canonical 
explicit symmetry violating terms, the canonical symmetry violating term of Eq. (2.3) is also renormalized. Because 
of the effects of explicit chiral symmetry breaking terms on this renormalization, the coefffcient mp of this symmetry 
violating term in the cutoff Hamiltonian need not be zero, even when the full relativistic theory has only spontaneous 
breaking of chiral symmetry. In fact, in the limit in which chiral symmetry is broken only spontaneously, the same 
constituent mass scale may appear in both the kinetic energy and the symmetry breaking interactions. This allows the 
quark constituent mass to set the scale for most hadron masses, and yet enables chiral symmetry breaking interactions 
to be sufficiently strong to make the pion massless. Thus light-front power counting criteria and renormalization allow 
us to introduce effective interactions which explicitly break chiral symmetry and yet may still ensure a massless pion 
in the relativistic limit g ^ g^. 

E. The Artificial Potential 

The task of solving the light-front Hamiltonian at the relativistic value of g is far too difffcult to attempt at the 
present time. Instead the goal of this paper is to define a plausible sequence of simpler computations that can build 
a knowledge base which enables studies of the full light-front Hamiltonian to be fruitful at some future date. 

A crucial step in simplifying the computation is the introduction of an artificial potential in the Hamiltonian. 

A primary rule for the artificial potential is that it should vanish at the relativistic limit. For example, the artificial 
potential might have an overall factor (1 — g"^ /g1) to ensure its vanishing dX g = gs- This rule leaves total fiexibility 
in the choice of the artificial potential since no relativistic physics is affected by it. 

The first role we propose for the artificial potential is that it ensure that the bound state structure of the theory at 
very small g is similar to the actual structure seen in nature. However, we also want the weak-coupling behavior of 
the theory to be similar to QED in weak coupling. This will ensure that methods of computation already developed 
for QED will be applicable. To ensure these connections we propose to structure the artificial potential so that bound 
state energies all scale as g"^ for small g, just as QED bound state energies scale as e^. We then demand a reasonable 
match to experiment when the scale factor 5^ is set to g^, even when higher order corrections (of order g^, 5*, etc.) 
are neglected. 

The second role of the artificial potential is to remove unfortunate consequences of the nonzero gluon mass from the 
gluon exchange potential. Due to the nonzero mass, the gluon exchange potential falls off too rapidly in the transverse 
direction, while being too strong in the longitudinal direction. In the longitudinal direction there is an instantaneous 
linear potential of order g^ which normally would be completely cancelled by the onc-gluon exchange; but due to the 
nonzero mass, the cancellation is incomplete. Without the artificial potential there are even instabilities that prevent 
the existence of stable bound states at weak coupling. More details on these instabilities are given in Section X.B. 

The third role we assign to the artificial potential is to give the weak-coupling theory a structure close to the CQM 
so that past experience with the quark model can be used to determine the precise form of this potential and to fit it 
to experimental data. To ensure this we will define an initial calculation which involves QCD complications only in 
a very minimal form. 

The final role of the artificial potential is to incorporate a linear potential in both the longitudinal and transverse 
directions to ensure quark confinement for any g. This is important for phenomenology. It is also needed for studying 



the roles of a linear potential and where it might originate, especially in the relativistic limit where the artificial 
potential vanishes. 

Only broad principles will be laid down here for the artificial potential. Its construction in detail will require a 
collaboration between specialists in the CQM and in QCD perturbation theory. 

The basic need is to incorporate the qualitative phenomenology of QCD bound states into the artificial potential. 
This qualitative phenomenology comes from three sources: kinetic energy, Coulomb-like potentials, and linear poten- 
tials. We propose that all three terms should be present in the weak-coupling Hamiltonian and all should have the 
same overall scaling behavior with g in bound state computations, namely g^. The kinetic and Coulomb-like terms 
can be constructed directly from the canonical Hamiltonian combined with a one-gluon exchange term, the latter 
obtained for the case of zero gluon mass. The Coulomb term, if represented in position space, has the usual form 
g'^/r, except that the definition of r we propose is 



{p'^Sx 



and there actually are two terms which are added: 



nin 



p+j, p'+r' 



5x\ (2.5) 



(2.6) 



where 5x~ is the light-front longitudinal separation of two constituents, 5x±_ the transverse separation, rric the con- 
stituent mass, and p+ the constituent longitudinal momentum. The prime refers to the second of the two constituents. 
The p^ /nic in the definition of r ensures that the dimensions match. The positive or negative SU{'i) charges must 
also be inserted. See Appendix C for details. 

A few comments regarding the form of the Coulo mb p otential are in order here. Our nonrelativistic limit is 5 — > 
and not the rric — > 00 limit studied earlier in Ref. ||35[| . In this limit p^ is held fixed while 5x~ scales as -j-. In 
light-front coordinates, the factor p^ is necessary for dimensions as already stated and in the nonrelativistic limit 
p^ is proportional to the center of mass momentum of the bound state independent of the relative coordinate 5x^ . 
Relativistically, r will need a careful definition to avoid possible disastrous behavior when p^ << p'^ or vice versa; 
this problem has not been studied. Also in the relativistic case, p'^ does not commute with 5x~ ^ so p'^Sx^ must be 
symmetrized to preserve Hermiticity. 

Finally, we need a linear potential term — terms proportional to r and r' . In Coulomb bound states, both are of 
order I/g^. Hence to achieve an energy of order g*, the linear potential must have a coefficient of order g^ . Thus the 
linear potential term, in position space, would be proportional to g^r. To be precise, and get dimensions straight, its 
form is 



V.-g^p\^ + ^l (2.7) 

with /3 a numerical constant. 

The linear potential needs to exist between all possible pairs of constituents: qq,qq,qq,qg,qg, and gg, where q,q 
and g stand for quark, antiquark and gluon respectively. The linear potential must always be positive (confining) 
rather than negative (destabilizing). We show in Section X.B. that it cannot involve products of SU{3) charges as 
the Coulomb term does. The Coulomb term could be given a Yukawa structure rather than the pure g^/r term. All 
of the potential has to be Fourier transformed to momentum space and then restricted to the allowed range of both 
longitudinal and transverse momentum after all cutoffs have been imposed. 

The artificial potential must also contain counterterms that remove the unwanted components of the one-gluon- 
exchange and instantaneous- gluon terms. Thus in addition to the order g^ linear potential, linear in both the longi- 
tudinal and transverse directions, there is a subtraction to remove the order g^ linear potential in the longitudinal 
direction. The potential removed is the potential that remains after the incomplete cancellation of the instantaneous 
potential in the canonical Hamiltonian by one-gluon-exchange terms. 

To ensure that the artificial potential vanishes at gs without destroying its weak-coupling features, the subtraction 
term has to be treated with care. We suggest the subtracted linear potential be multiphed by (I — g^/gf ) so that the 
subtraction begins to be negated only in order g^, which is smaller for small g than the artificial g^ linear potential 
that needs to be dominant. All other terms in the artificial potential can be multiplied by (f — g^ / g1) instead. 

To ensure that the Coulomb term shows Coulomb behavior, at least roughly, at typical bound state sizes, it is 
important that the mass used in any Yukawa-type modification of the Coulomb term scale as g^ rather than being a 
constant mass. 



Finally, the linear potential should be invariant to the full SU{6)w flavor-spin symmetry of light-front quarks and 
antiquaries discussed by Lipkin and Meshkov |^, leaving all SU{6)w violations to come from quark masses in the 
kinetic energy, Coulomb terms, and the more complex renormalization subtractions. This means that the quark mass 
scale appearing in the linear potential will be the same for up and down quarks. 

The artificial potential would be added to the cutoff canonical Hamiltonian. Then a complete set of counterterms 
would have to be added to ensure that the theory has a limit as the cutoffs are removed. 

The artificial potential is designed to allow confined few-body states to emerge from a field theory. The constraints 
that dictate its design are severe, and the type of potential we have discussed is the simplest we have found that meets 
these constraints and does not require large cancellations that thwart a constituent picture. However, this potential 
has serious flaws. Probably its worst flaw is that it confines everything including color singlets to a single region of 
space. Any confining interaction that is purely attractive between all particles does not produce scattering states. 
QCD must manage to produce strong attractive forces between color charges as they separate, without producing such 
strong forces between color singlets. In lattice gauge theory this is arranged through gauge invariance, which forces 
links to exist between color charges but not between color singlets. In a field theory calculation without manifest gauge 
invariance, there is no such obvious mechanism to turn strong forces on and off by hand. Higher order calculations 
will hopefully provide clues to how QCD produces the important phenomenological effects of our artificial potential, 
without producing the unwanted side-effects. 

III. LIGHT-FRONT POWER COUNTING: CANONICAL STRUCTURE 

To construct the bare cutoff LFQCD Hamiltonian with counterterms we will try to follow the procedure adopted 
in standard canonical covariant theory as much as possible. In the canonical theory, when we begin the analysis 
the bare cutoff Lagrangian is unknown since the counterterms are unknown prior to the analysis. One can perform 
perturbative calculations and determine the counterterms order by order in perturbation theory. If one starts with the 
canonical terms, which include all the possible terms in accordance with covariance and power counting, and if one 
maintains manifest covariance and gauge invariance at all stages of the calculation, one finds that the counterterms 
are also of the same form as the canonical terms. Thus in the standard procedure the only unknown parameters are 
just constants, namely, masses and coupling constants. 

Because we are interested in the cutoff light-front Hamiltonian, we cannot use covariance as a guide. Instead 
we have to consider power counting based on the limited kinematical symmetries of the light-front. As discussed in 
Section III.C, one consequence of the use of light-front coordinates is that counterterms may contain entire functions of 
momenta, which is of course a radical difference from renormalization in an equal-time formulation. Before discussing 
renormalization, however, we establish in this section the power counting rules for determining the possible structure 
of operators in the canonical Hamiltonian. 

A. Light-front power counting 

Light-front power counting is in terms of the longitudinal coordinate x~ and the transverse coordinate x±. Why 
are the two coordinates treated differently? The main objective of the power counting analysis is to deduce the most 
general structure of divergences that arise from increasing powers of the interaction Hamiltonian in perturbation 
theory. But power counting based on the kinematical symmetries of the light front is different from power counting 
based on the kinematical symmetries in equal-time coordinates. This is immediately transparent from the light-front 

dispersion relation for free particles, k~ = -*- _,_ . Because the energy factors into separate fc+ and k± dependencies, 

the subtractions are not constants. For example, when k± gets very large the energy diverges no matter what k~^ is. 
Thus, in general, we get a divergent constant multiplied by a function of k~^. In position space this translates into 
divergences at small x± being nonlocal in x~ and spread out over the light front. A similar result follows for the case 
when k~^ gets very small. This situation is to be contrasted with the equal-time case. Recall the relationship between 

energy and momentum in equal-time theory, E = \/ k^ + m? . In this equal-time form, if /c^ — > cx) while k^ is fixed, 
the kz dependence becomes negligible and arbitrary functions of k^ cannot arise. 

When we analyze the canonical light-front Hamiltonian, we confirm that indeed it does scale differently under 
x~- and x^-scaling (strictly speaking, a unique transverse scaling behavior holds only in the absence of masses). To 
determine the scaling properties of the canonical Hamiltonian it is useful to recall the canonical equal- a;+ commutation 
relations obeyed by the field variables in the two-component formulation p^,^- For the two-component gluon field 
variables, A^, i = 1, 2, we have 
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[Al{x),Aliy)]^,^^^ ^ -'-e{x- y-)5^x^ y±)S,,Sat. (3.1) 

For the two-component quark field variables ^(x) (the non- vanishing upper components of ^+(a;)) we have 

{6(x),Cl(2/)},+=,+ = S'ix - y)5,,. (3.2) 

Consider the scaling of the longitudinal coordinate 

x^ -^ x'^ — sx^ . (3.3) 



From the canonical commutation relations ( p.l| ) and (3.2), under this scale transformation A^{x ) -^ 
Ul{s)A^{x~)Ui{s) = A^{sx'~) and £,{x^) -^ U{ {s)£_{x^)Ui{s) ~ s^^^sx^). Ui{s) is the unitary longitudinal scal- 
ing operator. Next consider the scaling of the transverse coordinates 

x± -^ x'^ = tx±. (3.4) 

From (O) and (|3l), A\x^) -^ ul{t)A^{xi_)Ut{t) = iA^(te_L) and S,{xi_) -^ U}{t)^{xi_)Ut{t) = t^{txj_). Here 



Ut{t) is the unitary transverse scaling operator. Hence, the power assignments for the field variables are 



Ai^ : 


1 

Xl_ 


i ■■ 


1 1 


Vx- xi_ 


We also need the power assignments for the derivatives 




d± : 


1 

x± 


9+ : 


1 

x~ 


and the inverse longitudinal derivative 




1 

d+ 


x~ . 



(3.5) 



(3.6) 



(3.7) 



It can be verified that under longitudinal scaling, the canonical Hamiltonian density is invariant, Ti. ^ Ti' — Ti. 
(see Section IV), and the Hamiltonian 

H ^ H' ^ I dx'-(fx±H' = sH. (3.8) 



So the Hamiltonian H = P~ scales just like x", irrespective of whether there are masses present or not. Similarly, 
for the interaction Hamiltonian density in the absence of the helicity flip interaction, under transverse scaling H -^ 
H' ~ jsli. and thus 

H -> H' = 4-ff- (3.9) 

So in the absence of masses, the Hamiltonian scales inversely as x\. 

The Hamiltonian we will consider also includes mass terms. Under scale transformations masses scale as constants. 
The Hamiltonian density for the mass term in the free part, for example, scales as 7i — *■ 7i' = t?^- Thus the 
Hamiltonian does not have a unique scaling behavior when masses are present. For dimensional analysis we assign 

TO : — . (3.10) 



In the presence of masses, we assign 

n 



1 



H : \. (3.11) 



These assignments have meaning only in dimensional analysis. 
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B. Structure of the canonical Hamiltonian from power counting 

We assume that the canonical LFQCD Hamihonian density is a polynomial in the six components m, A^, S,, d^ , d±, 
and ^. Then the most general structure we can build for the canonical Hamiltonian density, which has dimension 

{^^^nA^,d^,m)^-^P{d+)-P. (3.12) 

Here the expression {A±,d±,m)'^~'^^ stands for monomials of order 4 — 2p in any combination of the three variables 
A±, d±, m. The resulting structure is 



P^O: Ai,Ald±,Aldl,m^Al 

1 
d 
1 \2 



P=^- ■7^{C^^){Al,A^d^,dl,mA^,md^,m^ 



9+ J (?^^)'- (3-13) 

We need not specify the action of d± since we assume that the inverse of d± does not occur in the canonical Hamil- 
tonian. But we need to specify the action of the integration operator ^^ , which occurs only because of the constraint 
equations that eliminate the dependent variables ■(/;_ and A~ . From -tp^ elimination we get the structure 

±{{^,^^){d^,A^,m)}, (3.14) 

and from A^ elimination we get the structures 

^{d^A^},{^na\A^d+A^}. (3.15) 

We will require that inverse powers of 9^ only appear in these specific combinations in the canonical Hamiltonian. 
Subject to this constraint, we can include both positive and negative powers of d^ . 

Specifying the precise way in which -J^ acts, we can enumerate those terms which obey our canonical rules: 

p = 0: {A^)\ {A^fd^, {A^)^d^f,m^{A^f, 

{A^d+A^)-^id^A^),iA^d+A^)i-^)\A^d+A^) 

P = 1: i^e)i^nA^d+A^), i^C^)±{d^A^), 

{m,d^,A^)^^-^{{m,d^,A^)^} 

P = 2: i^eK^rm^). (3.16) 

Comparison with the Hamiltonian derived by the standard canonical procedure (see Section IV) reveals that the free 
terms rn?A\ and m.C^9_L^^ are absent. The absence of the first term is because of the presumed gauge invariance. 
A term like the latter does appear in the free part when ^_ is eliminated, but it is cancelled by a similar term. When 
the smoke clears, no 7^ matrices appear in the free part of the canonical Hamiltonian. Note that our list excludes 
terms that are equivalent through an integration by parts. This integration by parts is justified because we assume 
cutoffs are present (see the discussion in Section V). 

IV. LIGHT-FRONT QCD HAMILTONIAN: CANONICAL TERMS 

In the previous section we motivated the structure of allowed terms in the canonical Hamiltonian based on power 
counting. In this section we use the standard canonical procedure to determine the explicit form of the Hamiltonian, 
including color factors. This gives a Hamiltonian with terms which obey the canonical rules established in Section 
III. These canonical rules allow us to extend our Hamiltonian to include a gluon mass term; and as we demonstrate in 
Section VIII, our cutoffs will demand such a mass term to appear through renormalization even if we did not include 
it in the starting Hamiltonian. We choose to extend the definition of the canonical cutoff LFQCD Hamiltonian to 
include a massive gluon from the outset. Of course, this is just the starting point, for renormalization will eventually 
force us to add to our Hamiltonian counterterms which do not obey the canonical rules. 
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A. Canonical Hamiltonian 

In this Section we follow Ref. M, to which the reader is referred for further details. The QCD Lagrangian is 



C = -^TiiPt^'^F^,) + ^Pih^D^ - mF)i,, 



(4.1) 



where F'"' = 9^ A" - d" A'' - ^g[A^', A''] are the gluon field strength tensors and A^' = Y.a ^a^" are the 3 x 3 gluon 
field matrices, with T" the Gell-Mann 5f/(3) matrices: [T", T''] = i/^fe^T^ and Tr(T'^r'') = ^Sab- The field variable 
ip describes quarks with three colors and Nf flavors, and D'^ — d^ — igA^ are the covariant derivatives, while mp is 
a Nf X Nf diagonal quark mass matrix. 

In light-front coordinates with the light-front gauge A^ ^ A^^ + A^ = 0, the Lagrangian can be rewritten as 



-C = \F:\d-A,) + #;(9-^+) --H- lA-Ca + \{^lC + CH-) 



(4.2) 



where 



n 



i {E-^ + B-^] + {vA{a± • (*5x + gA^) + Pmp}^_] 



-d+{E-A-)-d\ElA- 



(4.3) 



is the Hamiltonian density and 



Ca 



-d+E- - {d^El + gr'^AlEl) + gV^^rV, 



C — id^tp- — {ia± ■ d± + ga± ■ A± + (3mF)'4>+- 



(4.4) 
(4.5) 



Here, E- = -\d+ A' , 



S* = —^d^A]^ and Bp^ = F^^ are components of the light-front color electric field and 
the longitudinal component of the light-front color magnetic field, and ip^ and ip_ are the light-front up and down 
components of the quark field: ip = "0+ + V'-j V'i = ^^-0 = ^7^7^^, where A+ 4- A_ = /, A±^ = A± and A+A_ = 0. 
In (abelian or non-abelian) gauge theory, only two components — the transverse components — of the vector gauge 
potentials are physically independent degrees of freedom. From the equations of motion, it becomes clear that the 
independent dynamical degrees of freedom in LFQCD are the transverse gauge fields A^ and the up-component quark 
field ■(/'+• The Lagrangian equations of motion lead to Ca = and C = 0, which imply that the longitudinal gauge 
fields A^ and the down-component quark field ip- are Lagrange multipliers. Furthermore, in light-front coordinates 
the four-component fermion field can be reduced to a two-component field. The two-component quark field can be 
explicitly formulated in a light-front representation of the 7-matrices defined by |37|] 



7+ = 


■ 0' 

2i 


•) 


7 


= 


'0 



-2i' 



1 


f = 


-ia' 
ia' 


7 


7' = 


[a^ 1 
_ ~a\ 



(4.6) 



Then the projection operators K± become 



A^ 



1 




A_ 




1 



(4.7) 



and 



V'^ 



^_ 



{j^)[a'{id^+gA^) + imF]i 



(4.8) 



Hereafter, we shall simply let ^ represent the (two-component) light-front quark field. The canonical LFQCD Hamil- 
tonian becomes: 
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H 



_|_ y_ rabc tade j^i j^j j^i j^j 



S}{(Ji_ ■ {id± + gA±) - imp} 



X ( -^ ) {cr± • iidi_ + gAi_) + imp}^ 



gd'Al ( — ) ir'^^Aid+Ai + 2C^T'^0 



gW 1 



2 \d 



{/"'"'Ald+Al + 2^tT°^) 



If zero modes are retained there are additional surface terms |8| . 



(4.9) 



B. Counterterms to canonical terms 



Light-front infrared singularities arise when one eliminates the unphysical degrees of freedom by solving the con- 
straint equations, which are the source of the operator 1/9+ in (4.9). This operator produces tree level, light-front 
infrared divergences in the instantaneous four-fermion, two-fermion-two-gluon, and four-gluon interactions. These 
divergences require counterterms, which we now construct. 



1. Instantaneous gluon exchange 
First, we consider the instantaneous four-fermion term 



H. 



gggg 



-2g2 J dx-d'x^ I (e^T'^O (^) i^^T^O 



In terms of the quark color current j^"' = 2^ty°^ and its partial Fourier transform 

1 



we have 



J^^i^-,^!-) 



9 



2(27r) 



dp+jr(p+,a.^)e^P "" 



^....-g^y rf'^^y_^(p+)2 



rfT)+ - 

-^j+'^(p+,^±)jr(-p^'^^)- 



Thus we see explicitly that Hqqqq has potential divergences. 

In the cutoff Hamiltonian |p+| is restricted to be above e. That is. 



+ OC 



dp-^ 



dp^ 



dp"' 



To find the divergent part we can do a Taylor expansion of the integrand: 



ip 



The linearly divergent term is 



%''ip*.':' 



i^h°-^)-l^ 



p+=0 



(4.10) 



(4.11) 



(4.12) 



(4.13) 



(4.14) 
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l^^/dVjr(0'-^)3r 



dVJr(0,a;i)j+'^(0,a;i). (4.15) 



The logarithmically divergent term vanishes because of the symmetric cutoff. Let us suppose we had used separate 
cutoffs e+ and e_ for small positive and negative p+, respectively. Then the potentially logarithmically divergent term 
is 

'-^\n(^-±^ j Sx^ j dx- j dy-j+\x-,x^){x- -y-)j+%y-,x^), (4.16) 

which is finite if we choose e+ = /xe_, with /i some number. This particular operator vanishes because of symmetry 
under x^ —^ y~ exchange, but it forces us to reconsider our original neglect of zero modes in the solution of the 
constraint equations. 

We are finding a potentially divergent sensitivity of the Hamiltonian to modes with small longitudinal momentum, 
and we regard this as a signal that finite zero mode effects may arise. The potential divergence arises from the 
constraint equations and reflects contributions from A^ close to zero longitudinal momentum. It is extremely naive 
to suppose that the canonical constraint equation for A^ reproduces the nonperturbative effects of a zero mode in 
any simple manner. We discuss the divergences that arise from the exchange of small longitudinal momentum gluons 
with physical polarization below, but even here we can think of the counterterms produced by the constraint equation 
as arising from the exchange of small longitudinal momentum gluons with unphysical polarization. Once one realizes 
that an exchange is involved and that the divergence is independent of transverse coordinates, it becomes clear that 
even the assumption that these counterterms will be local in the transverse direction is naive. We therefore allow 
terms in the Hamiltonian of the form 

'-I- J d'x^dx- J dy-d'y^j+-{x-,x^){x- -y-)OG{x±-y^)j+^{y-,y^). (4.17) 

Og is a function of the transverse variables and is restricted by dimensional analysis, kinematical boost invariance, 
translational invariance, and invariance under rotations about the longitudinal axis. Oq must be odd under x±_ -^ y± 
to keep (4.17) from vanishing, and this is not possible here because gluon exchange gives only an even number of 
polarization vectors e\ with which to contract the transverse indices. Thus there is indeed no finite counterterm 
associated with logarithmic divergences in Hqqqq, and our discussion serves only to illustrate one way candidate 
vacuum interactions can be identified. 

As discussed previously, the requirement of boost invariance means that there can be no finite part associated 
with the counterterm for the linear divergence (4.15), whose form breaks longitudinal boost invariance. Thus the 
counterterm for the canonical instantaneous four-fermion interaction is 
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H^L = -^- / d'x^ / dx-j+\x-,x^) / dy-j+-{y-,x^). (4.18) 



In a similar way, the Hamiltonian counterterms for the instantaneous two-fcrmion-two-gluon and the instantaneous 
four-gluon interactions arc found to be 

Hqq'aa^ = -~-J <fx^ j dx-j+-{x-,x^) j dy-jj%y-,x^). (4.19) 

and 

H^q^^^-J d^^^j dx-jj^{x-,x^)j dy-j+^{y-,x^), (4.20) 

where the gluon color current is j+" = f"-^^A\d^A\. These three counterterms precisely remove the tree level light- 
front infrared divergences. 



2. Instantaneous fermion exchange 

There are several terms in the canonical Hamiltonian which contain a single inverse derivative ^. With one 
exception, all these terms can be rewritten so that ^ acts on a single field operator. Such terms do not give rise 
to divergences since the zero modes are removed from single operators. The exception is the two-quark-two-gluon 
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interaction involving an instantaneous fermion exchange Hqqggi . Here ^ acts on a product of field operators and 
hence can give rise to a logarithmic divergence (see also Section V). 
Denoting (7± ■ A±[x^ , x±) S,{x^ , xj_) by f{x^, x±), we have 

Hqqggi = 5^ / dx^(fx±p{x^,x±)-—^f{x^,x±). (4-21) 

Introducing the partial Fourier transform as before 

/oo 
dx~eiP^-'f{x-,x^), (4.22) 

-OO 

we find 

H,,ssi = f fd'x^ r '^f\p+,x^)f{-p+,x^). (4.23) 

A potential logarithmic divergence disappears if we use a symmetric cutoff. 

Once again, though, we must be careful with this logarithmic divergence. While there is no divergent dependence 
on the symmetric longitudinal infrared cutoff e since the divergence due to small p^ is cancelled by the divergence 
due to small — p'*', we still have to worry about non- vanishing contributions from states with exactly p'^ = — that 
is, from exchange of zero-mode fermions. Once an infrared cutoff is employed, we have eliminated the possibility of 
computing this contribution; and so we should include in the Hamiltonian terms which can counter the effects of this 
exclusion. Such counterterms have the form 

H^q'rjgi=9^ j dx-d^x^ j dy-d^yU\x-,x^)OF{x^-Vi_)J{y-,Vi_), (4.24) 

where dimensional analysis reveals that Op scales as 1. Since this term comes from zero- mode fermion exchange. Op 
is a function of transverse variables and the Pauli matrices a^. Since this term arises from vacuum effects, it is not 
restricted to have a polynomial dependence on transverse variables and the rcnormalized constituent mass. Thus a 
term like Op ~ m~j^ a± ■ d± is allowed and should be included in the Hamiltonian. Such a term explicitly violates 
chiral symmetry and need not be small in the relativistic limit. 

This is then the first point at which an arbitrary function enters the Hamiltonian. At the simplest level, the 
function Op(xj_, yj_) must be determined phenomenologically by fitting bound state properties. It is an open question 
whether theoretical techniques such as coupling coherence will determine it a priori, whether the requirement of 
relativistic invariance will determine it completely a posteriori, or whether phenomenology will still be necessary 
when renormalization is done at higher orders. 

C. Canonical Hamiltonian: free plus interaction terms 

Now we add explicitly the gluon mass term to the LFQCD canonical Hamiltonian together with the counterterms 
for instantaneous interactions. The Hamiltonian is written as a free term plus interactions 

H^ f dx-d^x^{Ho+n^nt), (4.25) 

with 

Ho = l(d^Ai){d^Ai) + \mlAiAi+e ( '^^ + '"^ ) ^ 

li-int ^ Va + Hqqg + iLggg + llqqgg + llqqqq + li-gggg. (4.2d) 

Here Va is the artificial potential, and 

+ cr_L • A_L f T^ j ((T± • 9_L + mp) 
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n^sg = gr"' S^d'AiAlAi + (a^l) (^^) {Ald+Ar)'j 



(4.27) 



(4.28) 



'^qqgg = 9^ i ^^o"i • ^i ( T^ ) o-± • A^^ 



1 

i9+ 



2gU{r'>^Ald+Al) (^^y{i^T'^0]+Kli,i+HgL2 



^ "qgggi ' "qqgg'^ 



n,,,, = -2g4(etr'^e) (^)'(Ctr''o| +i?f,I, 



(4.29) 
(4.30) 



T-/ — d^fabcfadel A^ AJ A^ A3 

'^gggg ^ . J J \ ^b^c^d^e 



-2{Ald+Al)[ — ] {A^^d+Ai)}+H: 



CT 

gggg 



^ '^gqqql + 'I, 



99391 ^ rLgggg2- 



(4.31) 



In a careful analysis all of these couplings must be allowed to renormalize separately and the quark masses in Tig 
and Ti-qqg are not the same. For simplicity we use one coupling and one mass in our examples. It can be seen that 
all canonical terms are predicted from the light-front power counting (3.16). The artificial potential has the form 
described in Section II. E. The Coulomb part will involve the quark and gluon color currents j'l^" and j,j"" defined 
above, and the linear part will involve the color singlet currents j+ = 2^^^ and j'+ = AJj9+A^. 

As we have emphasized before, the division of the bare Hamiltonian into a free part and an interaction part is 
arbitrary. We have chosen the free part to be that of massive quarks and gluons since self mass corrections are needed 
in any case. Note that these terms are allowed by the power counting criteria. Our computations will involve two 
stages. The first stage is the computation of the effective Hamiltonian, which is to be performed perturbatively. The 
artificial potential will be considered as part of Tii for this stage. The second stage is the determination of bound 
states, and as discussed before we will want to include the (nonrelativistic) Coulomb and linear parts of Va as part 
of the unperturbed Hamiltonian for this stage. 

V. POWER COUNTING: PHASE SPACE CELL ANALYSIS AND THE STRUCTURE OF 

DIVERGENCES 

A. Phase space cell division and field operators 



Here we perform the power counting analysis of ultraviolet divergences and infrared divergences arising from the 
products of interaction Hamiltonians. To facilitate this analysis, we first introduce a phase space cell representation 
of field variables in terms of wave packet functions [ p8p9| ] . A realization of these wave packet functions is given by 
wavelets. The present considerations are only qualitative in nature, and hence we do not need their many interesting 
mathematical properties pO| ; however, they may be needed by anyone attempting the nonperturbative diagonalization 
of the renormalized Hamiltonian. 

The motivation for a phase space cell analysis and an example of phase space cell division is provided in Appendix 
D in the context of the more familiar case of equal-time Hamiltonians. Here we discuss the division of phase space 
cells in the light-front case. For the "plus" momenta we choose 



Vtiax >Pt >P2 >P3 



^ r^min 



(5.1) 



2J- 



with p^,i — ^Pi ■ For the transverse momenta we choose p±Q = 0, p±i = m, and then a sequence pj 
until the cutoff p±max is reached. Note that we always use a logarithmic scale for momenta. 

For a momentum space cell with centroids pj and p±j which is an annulus in the transverse momentum, position 
space is divided into a linear grid of cells all with width Sx^ — -^ and Sx± = :::^. These cells are labelled with 



1 

p7 
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position indices I and l± = (P, l^). A phase space cell is denoted by the indices ijl = {i,j, I , P, Z^). A wave packet 
function that belongs to the complete, orthonornial basis set is denoted in position space by <j)iji{x); and its Fourier 
transform is denoted by (piji {p) . The x-space widths of these functions arc 6x^ and Sx±j and they are assumed to 

vanish outside of the cell. Normalization implies that (t>iji is of order /,."'" , while 4>iji{p) is of order \ 5x^ 6x±j. 

Next we consider field operators. For the quark field operator, we have 



i{^) = Y.'^ 



dk+d^kA 

2(27r)3 



5fe,,e-*^+4.,e*^ 



with 



Thus 



{bk,s,bl,^,^,} = {dfc,„4,^^,} = 2{2TTfS{k+ - k'+)6^{k^ - k'^)6s 



6, d 



' X Xj_. 



(5.2) 



(5.3) 



(5.4) 



We wish to carry out the analysis in such a way that the creation and destruction operators are of order unity. To 
achieve this goal, we introduce phase space cell operators 



and 



where 



Thus 



lij 



ijl 



bijis , d^jis ^ 1 and 4iiji{k) ^ J x^ xj_j 



6 (a;) = J2 [b^jis't'^jii^) + 4jis<f'Ui'> 

ijl 

It follows that the order of magnitude of ijl*^ mode's contribution to the operator ^ is 
Next consider the transverse gluon field operator. We write 

dk+d^k. 



/SX- ^^±3 ' 



AAx)^J2 



2Vfc+(27r)3 



^ r, ^—ik-x J. * t „ifc-: 



with 



K,A,4',A'] - 2{2TT)^6{k+ - k'+)5\k^ - k'^)S^x'. 
The operators ak and aj, are again unsuitable for our purposes. We define 



ak,\ = /^4>iji{k)aiji\ 

ijl 



where (j)iji{k) ^ y ^^i Sx±j. Now, 



, , , V- /■ dk+d^k^ 



e±A^jj7(fc)ay/Ae '''■^ + e\)^(i)lu{k)a\^ixe'^' 



(5.5) 



(5.6) 



(5.7) 



(5.8) 



(5.9) 



(5.10) 



(5.11) 



(5.12) 
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The order of magnitude of ijl* mode's contribution to A±{x) is 337— • 

Because of the cutoffs we impose, the wave packet functions (fiiji (x) do not have zero modes; that is, / dx(j)iji {x) = 0. 

We must also specify the order of magnitude of 9^, 9+, and (9+)^^ when acting on a specified subset of wavelets 
in an operator product. For qualitative purposes each gradient needs to be applied only to the most rapidly changing 
wavelet in the subset, so that d±_ is of order {5x\_)~^ , where 6x^ is the smallest width in the subset. A similar result 
holds for d^ . The effect of the inverse operator {d^)^^ — which involves an integration — is also determined by the 
smallest width 5x~ in the product. Qualitatively, it amounts to multiplication of the subset by the width 5x^ , for 
the other, broader widths act as constants over this width. However, there is a special case we must note. 

To see this, consider the action of (9"'")^^ on a product of field operators. Now the wave packet functions associated 
with destruction operators have positive A:+ only, and the wave packet functions associated with creation operators 
have negative k'^ only. If the product contains only creation operators or only destruction operators, then the operation 
of {d^)^^ qualitatively amounts to multiplication by the smaller width among the product in accord with the rule 
above. But if the product contains at least one creation and one destruction operator, then wavelets with Fourier 
modes of positive k^ and wavelets with Fourier modes of negative k'^ will both occur in the product. When they 
have the same mean value of fc"*", the product acquires a zero mode and then (9+)^^ no longer produces a localized 
function but rather one which covers the whole space. Thus the action of a single (^) can lead to a logarithmic 
infrared divergence and the action of (gr)^ can lead to both linear and logarithmic infrared divergences already in 
the canonical terms. We have discussed these divergences and the counterterms which remove them in Section IV. B. 

B. Ultraviolet divergences and counterterms 

Consider ultraviolet divergences arising from the products of interaction Hamiltonians. Each term in the interaction 
Hamiltonian involves an integral over the product of three or four field operators. Each field operator is expanded in 
terms of wave packet functions. First, we consider the ultraviolet divergences arising from a product of two interaction 
Hamiltonians. We shall consider the special case of just two wave packet sectors, one being fixed and the other being 
arbitrarily large. This corresponds to a single divergent loop integral in the familiar language. The fixed wave packet 
width corresponds to the fixed external momenta in the loop diagram, and the varying wave packet width corresponds 
to the diverging loop momenta. 

We denote the widths of the fixed and varying wave packet sectors by {5x~ , 5x±^ and {5y^ , Syj_), respectively. 

The energy of the varying wave packet sector scales as A-- The energy diverges as Sy± -^ 0, corresponding to 

the traditional type of short-distance ultraviolet divergence, which we consider in this section. The divergences as 
Sy^ -^ 00 are light-front infrared divergences, which we consider in the next section. Consider ultraviolet divergences 
which occur when Syj_ <C Sx±. We set Sx~ = Sy~ . The interaction Hamiltonian contains only cubic and quartic 
terms in the field operators £.\ £,, and A^. We now observe that, as a consequence of momentum conservation at the 
vertices, to produce an ultraviolet divergence at least two of the three or four operators have to refer to a y-sector. 
We also note that a transverse gradient will yield jJ— unless it applies to a product of x-sectors only. 

Now consider the terms in the interaction Hamiltonian. If all operators refer to y-sectors, each term scales as tt- 

except for the mass term (the helicity-flip interaction), which scales as -f^- But if some operators refer to x-sectors, 
there will be additional factors of -^^. For each field operator (^^, S,, Aj_) belonging to an x-sector and for each 
derivative 9x acting on only x-sectors, there is a factor of -g^. This gives the contribution from a single sector with 
width 6y~,6y± as dyj_ gets small. But many wavelets with different Syj_ can contribute, namely, all wavelets with 
different centers that overlap a given x- wavelet (see Fig. 13). Since the transverse space has dimension two, the 
number of independent wavelets of linear size Sy± that overlap with the wavelet of size Sxj_ is of order (Sxj_/5y±)'^. 
Next consider the second-order perturbation theory formula 

Hi^Hj, (5.13) 

where E.D denotes the energy denominator. For qualitative purposes, we can replace the energy denominator by the 
energy of the y-sector, namely, tt-- Thus the contribution in second order of a given term to the effective Hamiltonian 
resulting from wave packet functions of a given width 5y±^ scales as 

(Ii)^(|^)^(lr)^(|r)^^(l7)" <-»' 
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where n counts the number of ^^ , ^, A±, and d± factors that refer to x-sectors. If n is greater than 4 there is no 
divergence, while for n less than or equal to four a divergence arises. When n is equal to 4, the contribution from 
the sector of a given size Sy± is a constant. Now we have to sum over all scale sizes Syj^. The scale sizes 6y± allowed 
are inverse powers of 2 all the way down to a minimum size -^ . A sum over all scale sizes down to -^ yields a term of 
order {InA /ln2), which gives a logarithmic divergence as A ^ oo. 

Thus the power counting rule for ultraviolet countertcrms is the same (in Sxj_) as for the canonical Hamiltonian; 
namely, that there can be at most four operators scaling as SxJ^ . Even though we did the analysis for two -ff/'s, it 
holds for any number of i?/'s, where the first Hj creates one or more y-type constituents, intermediate _ff/'s maintain 
at least one y-type constituent, and the last Hj destroys all the remaining y-type constituents. (We arc ignoring 
contributions "disconnected in y" — that is, products which contain an Hj with only x-sector constituents. These 
diagrams have already acquired countertcrms.) Each energy denominator may be replaced by the energy of the y- 

sector -7^ since each intermediate state contains at least one y-type constituent by assumption. If there are m factors 

Hi, there arc m — 1 factors j^. Each Hj scales like -^^ times powers of {Sy±/Sxj_). An overall factor {6xj_/Syj_)'^, 

counts the number of y-type wave packet functions that overlap with a single a;-type wave packet function. Thus the 
final rule for the scaling versus Sx± and Sy± is identical for m i//'s as for 2 _ff/'s — namely, divergences occur for 
products in the effective Hamiltonian of a:;-operators up to fourth order in ^^, ^, A±, and 9^. 

C. Infrared divergence and counterterms 

Light-front infrared divergences are large-energy divergences as can be seen from the light-front dispersion relation 

2 I 2 

p^ _ p± m ^ ^}jjq}j blows up as p+ -^ 0. We have already discussed the counterterms for light-front infrared 
divergences that arise in the instantaneous interactions in the Hamiltonian. In this section we discuss the infrared 
divergences which arise in products of the interaction Hamiltonian. 

As in the analysis of ultraviolet divergences we will separate out the divergences caused by a specific y-type wave 
packet function from the operator structure associated with an x-type wavelet. Let us discuss the differences from 
the analysis of ultraviolet divergences. First of all, the y-wavelet energy scale -^ is blowing up because Sy~ is 

y± 

approaching oo. Thus it will be assumed that 5y^ 3> 6x^ , while to avoid confusion with ultraviolet problems we 
assume Sy± is of order 6x±. In contrast to the ultraviolet case, light-front infrared divergences are caused by a single 
y-type wavelet with a width Sy^ . The reason for this is that only one y-type wavelet overlaps with a given x-type 
wavelet, since 6y~ is much wider than Sx~ . In this case momentum conservation at the vertices implies that to 
produce a light-front infrared divergence at least two operators in Hj have to be x-type. Each term in Hj scales as 

TT" times a factor: a) (jf^)^ for each ^ or ^^ referring to an .x-type, b) (^f^) for every d^ acting on an x-type 

' Sx 

fS 



operator and a reciprocal factor (|^) for every (9+) ^ operator when applied to a product including at least one 
x-type operator, and c) (f^) arising from the range of integration in each Hj involving an x-type operator. As in 



the case of ultraviolet divergences, in the effective Hamiltonian a term with no factor (^^) diverges logarithmically 
due to a sum over different widths 6y^ . However, exceptions to the above rules occur for terms in the canonical 
Hamiltonian involving (g+ ) or (q+)^ operators, which we will note below. 

Let us now catalogue the factors associated with the scaling behavior of various terms in the interaction Hamiltonian. 
If we ignore 9+ or -^ factors, then the rules are simple: an x-type ^ produces a factor (ff^)^^^, an x-type A±_ produces 
no factor, and the presence of any x-type operator in Hj produces (f^) as result of the change of integration width. 

Thus, ignoring 9+ and ^ factors, we can compute the impact of any product of x-type operators on the scaling 
behavior of Hj : 



cc 


no factor 


^±e 


{6x-/Sy-)i 


^±ee 


no factor 


A±A^ 


iSx-/dy-) 


A±A^^ 


(Sx^ /Sy^)2 


A^A^CC 


no factor 


^?? 


iSy-/Sx-)-^ 
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^^^^ 



{Sy /Sx ). 



(5.15) 



Next we classify again the possible x-type operators in Hj by the degree of divergence in Hj that results, when the 
effects of 9+ and ^ factors are included. Here it is assumed that d~^ or (1/9"'") is applied to a product that contains 
at least one x-typc operator. It is useful to recall the structure of terms in the canonical Hamiltonian (see Section 
IV), which we consider one by one in the following. 

When all the operators in the Hamiltonian refer to the y-sector, the Hamiltonian scales as Sy~ . We have to 
determine the scale factor that arises as a result of replacing two or more y-type operators by x-type operators. Let us 
first consider some specific examples. In the following we denote the operators that belong to the x-type by a subscript 
(x) and those belonging to the y-type by a subscript (y) respectively. Now consider the term in the Hamiltonian 
J dz~(Pz±^'^ A±-^^. Here we write the integration variable as z in order to distinguish dz from the widths 6x and 
6y. 



1 . For the case / dz d^z^(Aj^)(j,)^Jx^^(j.), from our rules we get (Jy ^ \j^ 

(j^) (arising from the range of integration). Thus the scaling behavior is like 5x^ which we write as 5y^ x 



X (|al)(fortwoe's) X (|f^)(forgV 



(|^). Hence the scale factor is (|^). 



2. For the case J dz^d'^z±{A±S_'^)(jj.)-^^(y) we get from our rules Sy^ x 
the range of integration) . Thus the scaling behavior is like (6x^ 6y 
Hence the scale factor is (j^)^ ■ 



(ll^)' (for^''') X (|f^) (arising from 
")2 which we write as Sy~ x (|^)2 . 



Now there may be exceptions to these rules when (d^)^^ or (9+)^^ acts on a product of operators, for then it would 
pick out a zero mode as described above. However, we have eliminated this possibility by imposing a longitudinal 
cutoff e on such terms and then adding a counterterm to remove the e-dependence. Because of this, integration by 
parts is allowed and can be used to determine the order of magnitude of (9+)^^ or (9+)^^. The exception to the 
usual rule occurs, then, when there are not x-sector wavelets on each side of the inverse gradient, for integration by 
parts allows us to re-write the product so that the inverse operator acts entirely on y-type wavelets. Then it picks 
out a y-sector width and so does not produce the factor {Sx~ /Sy~) or {dx~ /Sy^)'^ expected from the above rules. 
Thus, for example, consider the term J dz~d'^z±{^'' A±)ty\-Q^{^A±)(^xy With the zero mode eliminated, we are free 

to interchange the operation of g^ to get J dz^d'^z±-^(^{^^ Aj_)(^y-)){A±^)(^x)- Thus the scaling is 6y^ x (^f^)^ 
(from ^) X (t^) (arising from the range of integration). Hence the scaling behavior is {dx^ 5y^)^ ^ which we write 



(If-)^; 



and so the scale factor is (1^) 



as 5y X 

In the following we determine the scale factor for various terms in the interaction Hamiltonian when two or more 
operators in them are replaced by x-type operators, the remaining operators being of y-typc. Each -g^ may apply 
to y operators too, e.g, ^^ might mean -g+{A^^) with A^ of y-typc. The scale factor is unchanged by any such 
embedded y operators. 

Products of ^ and S, : 



c(a+)-^c 



no factor 

no factor (exception applies) 

{5x^ jdy^) 
{5x-/5y-f. 



(5.16) 



Products of A±_ and ^ : 



AU ■■ 


(Sx^/Sy-)^ 


id+)-\A^O 


(Sx^/Sy^)^ 


ad+y\d+A^) : 


{5x^ /6y^)2 


{d+A^id+r^c ■■ 


{5x-/6y-)i 


ao+r^A^ : 


{6x-/5y-)i 


A^{d+r'i : 


{5x- /5y-)i 


ad+y'A^ : 


{5x- /5y-)i 



(exception applies) 



(5.17) 
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Products of ^, ^, and A± 



Products of A± and A± 



e{d+)-^d+A^ 

A^{d+)-'e 

{d+A^)id+)-^e 

eid+r'A^ 



A^Ai_ 

A^id+)-^Ai_ 

A^id+)-^d+A^_ 

{d+A^){d+)-^d+A^ 

(d+r^A^d+A^) 

A^d+A^ 



{5x^ / 5y^ 
{Sx^ /Sjf 
(5x^ /Sjf 
{6x^/6y' 
{6x~/6y~ 
{6x^/Sy' 
{6x^/Sy' 



[Sx /Sy ) 

{6x-/Sy-f 

{6x-/5y-f 

{5x^ /5y^) 

no factor (exception applies) 

no factor 

{6x-/5y-f. 



(5.18) 



(5.19) 



Products of S,^ A^, and A±^ 



^{d+)-^A^d+A^ 
A^^d+r^A, 
A^{d+)-^iA, 



{Sx^/Sy' 
{6x^/Sy' 
{6x^/dy' 



(5.20) 



Products of ^, ^, and S, 



e(a+ 

e{d- 



r'e 



{Sx /6y' 
{6x^/dy' 



(5.21) 



The above list splits into four types: a) Products which contain no |^ factors. They are all components of the color 
charge density, either the quark or the gluon component. These products can occur in any number of iJ/'s without 
reducing the divergence of the overall product. If these products are the only source of x-type operators, the overall 
divergence behaves as Sy~ , that is, linearly, b) Products which contain a factor of {6x~ /Sy^)^ . These products 
involve a single ^. Since only an even number of ^'s can appear in the effective Hamiltonian, two such products 
are required from two separate i?/'s; and the result is to generate a factor {5x^ /Sy^), which implies a logarithmic 
divergence. If more than two such products occur there is no longer any divergence, c) Products which generate a 
factor {5x~ / 5y~). If any such product occurs once, a linear divergence is converted to a logarithmic divergence; and 
if they occur more than once, there is no divergence. Finally, d) products that scale as a higher than linear power of 
{5x~ / 5y^). Such products cancel any divergence — even if present linearly. 

As we have already explained, counterterms for linear divergences cannot have finite parts. Thus we are interested 
in logarithmic divergences. Even in this class, the power-counting analysis may be an overestimate since, for example, 
as seen in Section III.B, a logarithmic divergence may disappear when a symmetric cutoff is used. 

D. Examples 

Now we provide examples of the determination of the divergence structure in perturbation theory using the quali- 
tative phase space cell analysis described above. We consider the second-order shift in the energy of a gluon coming 
from a two-gluon intermediate state. 

As candidates for iJ/, we take 
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^/(i) = 9.r' J dz-d'z^d'AiAlAi, (5.22) 

Hm = -gr"' J dz-d'z^[j^d^A:)Ald+Ai. (5.23) 

We separate the field operators into low- momentum parts which contain wavelets of width Sxj^^Sx^ and high- 
momentum parts which contain wavelets of width Sy±, Sy^ . In the discussion below we drop factors of g and color 
factors. 

1. Ultraviolet divergence 

To avoid confusion with infrared divergences, we set 6x^ — 6y^ . For a candidate Hi, we choose -ff/ri). Remember 
that to produce an ultraviolet divergence at least two operators have to belong to the high sector. 

a) Consider -^/(i) ~ J dz~ d"^ z±{d^ A^ A^)(y^A-ls and the energy shift AE « H ^i^ ^^-j^ H j (^ly Referring back to the 

scaling formula (5.14), n = 2. Hence the scaling behavior of the energy shift is vj-^!-^ x (Sy^)"^ x ,g ^ .-i x {j^^)"^ x 
(sx^)^' ^^^^ is, AE w jj-^-yi, a quadratic ultraviolet divergence. 

b) Consider -ff/(i) ~ J dz~d'^z±d'^A-l.{A'^A^)(^y) and the energy shift AE « _ff/(i)— !-^i//(i). In this case n = 4 and 

hence the scaling behavior of the energy shift is p^ x iSyj_f x -^j^ x (ff^)^ x {^)^, that is, AE ^ ^^^^^^ , 
a logarithmic ultraviolet divergence. 



2. Infrared divergence 

To avoid confusion we set Sx± = Sy±. Remember that to produce an infrared divergence at least two operators 
have to belong to the low sector. Thus, for example, we may take 

J^/(i) ^ J dz-d'z^d'Al^{A'A')^,y (5.24) 

F,(2)« /'dz-d2zi(^a^^))(A^a+A^")(,). (5.25) 



a) Consider the energy shift AE w -H^/(2)"ET)-^-f(2)- According to our rules, -ff/(2) scales like Sy^ . The energy 
denominator produces a factor j^. Another -ff/(2) produces a factor dy^ . Thus AE w Sy^ , which results in a linear 
infrared divergence. 

b) Consider the energy shift AE w i//(2)-g^-ff/(i). According to our rules, ff/(i) scales like Sy^ x (^^). The 
energy denominator produces a factor j^. .ff/(2) produces a factor Sy^ . Thus AE w Sx^ , which results in a 
logarithmic infrared divergence. 

To complete the analysis one must determine what happens when Sy± <ti Sx± and simultaneously Sy^ 3> 6x^ , 
as well as other double orderings involving drastically different transverse and longitudinal widths. We have not 
completed analyses for all such cases. 

E. The structure of counterterms: summary 

We have found that the power counting rule for ultraviolet counterterms is the same (in 6x±) as for the canonical 
Hamiltonian, that there can be at most four operators scaling as j^ . However, the counterterm has a complex 
nonlocality in the longitudinal coordinates. The reason for this is that the longitudinal distance scales for x-type and 
y-type wavelets overlap even when ultraviolet transverse divergences are produced. Thus the ultraviolet counterterms 
are built out of products up to fourth order in ^, ^', A±, and d± (see Eq. (3.16)) and have completely arbitrary 
longitudinal structure. As a result of the constraint from longitudinal boost invariance only vertex counterterms can 
involve a priori unknown functions of longitudinal variables. 
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The counterterms for infrared divergences, on the other hand, involve arbitrary numbers of quark and gluon 
operators. They have a complex nonlocality in the transverse variables. This is in contrast to the divergent counterterm 
for the canonical instantaneous four-fermion interaction which is local in the transverse direction. If we identify the 
counterterms arising from infrared gluons (small longitudinal momentum) as the source of transverse confinement, 
then the unknown nonlocal transverse behavior would have to include confining effects at large transverse separation. 

From a physical point of view the appearance of long-range many-body interactions involving arbitrary numbers 
of quarks and gluons is inevitable if one allows a confining two-body interaction. It has long been appreciated in the 
CQM that if one uses confining two-body interactions, unphysical long-range van der Waals forces inevitably arise 
above 1-1-1 dimensions ma. In the CQM various schemes allow the potential to distinguish between colored objects 
in different hadrons, but we have no such operators at our disposal in a field theory. Therefore the only way to cancel 
the residual long-range multipole interactions between color singlets is with long-range many-body interactions. It 
would be a failure of our approach if such operators were not allowed by the same arguments that lead to a confining 
two-body interaction. The possible impact such operators have on the predictive power of the theory is not known 
yet. 

There is a question of whether the artificial potential will drastically alter the analysis of infrared divergences. We 
do not study this problem here. 

As mentioned in Section II, one might worry that the appearance of functions of momenta in the counterterms could 
destroy the predictive power of the theory and lead to non-renormalizability. However, these functions are needed to 
restore Lorentz covariance and gauge invariance as g —^ gg, and without such counterterms physical quantities will 
not even approach finite limits as cutoffs are removed for any coupling. 

We do not expect any new parameters to appear in the properly renormalized theory; so all new counterterms must 
actually be determined by the finite set of canonical free parameters. This problem has been studied in Ref. [ p2| , 
where a set of relationships called "coupling coherence" [Q were used to fix this relationship. Coupling coherence 
rests on the observation that only canonical variables should run independently with the cutoffs in renormalization 
group equations. To lowest orders in perturbation theory it has been explicitly demonstrated that coupling coherence 
fixes the bare Hamiltonian in scalar field theory, and that Lorentz covariance is restored by the resultant counterterms 
[ p9[ . Whether coupling coherence helps solve the infrared problem has yet to be demonstrated, but it is at least able 
to deal with spontaneous symmetry breaking [ p2[ . 

VI. A CUTOFF SCHEME FOR THE HAMILTONIAN 

In the previous sections we presented a qualitative analysis of divergences based on wavelets. Now we set up 
a precise momentum space framework for renormalization of the LFQCD Hamiltonian. We employ the "similarity 
renormalization scheme" |l4,^q] , which allows one to avoid the difficulties often encountered in traditional perturbation 
renormalization schemes when using plane wave states to study the structures of divergences. Thus, although we use 
plane wave states to determine the form of the effective Hamiltonian, we may still set up the renormalization to 
follow the qualitative wavelet analysis of Section V. Of course, the assumption here is that the artificial Coulomb and 
linear potentials will not affect the structure of the divergences. As we have seen, the canonical Hamiltonian already 
leads to divergences at the tree level due to k~^ getting small. Beyond tree level, products of Hamiltonians also lead 
to divergences both due to fc"*" getting small and k± getting large. We need to introduce a regulator before we can 
investigate these divergences and construct the corresponding counterterms. This is done here. 

A. Cutoffs on constituents 

Since the final step of the renormalization process involves the nonperturbative diagonalization of the Hamiltonian, 
we need to develop a cutoff procedure that is applicable to the Hamiltonian as a whole instead of just to perturbative 
calculations order by order. One way to accomplish this is to cut off the single particle momenta appearing in the field 
variables themselves. Cutoffs on the constituent momenta k^ and kj_, which we call "constituent cutoffs", obviously 
violate the longitudinal and transverse boost symmetries of the light front. What docs this violation mean? The 
direct effect of the cutoff is to eliminate states that would be present without the cutoff. For a given center of mass 
momentum, a key parameter is the lowest mass state that is not eliminated by the cutoff. This lowest mass depends 
on the center of mass momentum chosen, and is degraded if P± increases or P'^ decreases. 

If constituent cutoffs destroy light-front boost invariance except for a limited range of center of mass momenta, 
then, why do we want to use such cutoffs? One could instead employ — as has been suggested elsewhere, see e.g.. 



9Lp5 43 44] — cutoffs which preserve both transverse and longitudinal boost invariance. We shall call such cutoffs 
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"Jacobi cutoffs" for they act on the Jacobi or internal momenta of a constituent, which are defined relative to the 
center of mass momenta. Unfortunately, Jacobi cutoffs inevitably refer to extensive quantities of a state and thereby 
introduce nonlocalities in the Hamiltonian. If the cutoff on constituent i depends on the total momentum of the state, 
then a dependence on all spectators j will be introduced in any matrix element. Thus any counterterm may change 
when another spectator is included, which greatly complicates the renormalization process p9| . So there is a stiff 
price to pay for insisting on a regulator which preserves explicit boost symmetry. 

We prefer to employ a constituent cutoff scheme because of the conceptual simplicity of its implementation. This 
is especially important for LFQCD, where the final diagonalization of the Hamiltonian will be done on a computer. 
In addition, we will be encountering nonlocal effective potentials which result from the similarity renormalization 
scheme, and it will be important to avoid nonlocalities from the regulator itself. Choosing constituent cutoffs ensures 
there is a limited possibility of confusing between different sources of nonlocalities. Moreover, as we show below, one 
can choose constituent cutoffs which would allow boost invariance to be maintained to a good approximation within 
a large domain of center of mass momenta. In the next section, we discuss the details of the cutoff scheme and map 
out the center of mass domain (the values of P+ and P^) in which when a cutoff constituent appears in a state, the 
internal mass of that state is guaranteed to be above an effective cutoff A^. 

B. Details of the cutoff scheme 

We start with the formula for the total mass of an eigenstate of Hq (where Hq is the free part of the Hamiltonian) 
of n constituents with momenta [k^ , ki±) and masses rrii: 



A'P 



^^E 






(6.1) 



where P+ = ^^ k^' and P± = ^,j ki±. If we introduce the internal or Jacobi momenta Xi — -^ and gix = kn_ — XiP±, 
then 



m'=j: 



5^ 



(6.2) 



Thus the square of the total mass of the system is simply a sum of internal mass terms, which is a manifestation of 
the kinematical boost symmetries. 

We need to prove that the partial sums of internal mass terms can always be lowered by replacing two constituents 
by a single constituent with the same total q± and x and the lower of the constituent masses; that is, 
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(6.4) 



But the minimum value of ^ + ^ is {mi + 7712)^- Hence 
mi + qlj_ 



mj + g|_^ ^ {nil +m2y + (gi_L + q2±) > ^ 



Xl 



X2 



Xl + X2 



+ {qii- + 92_l)^ 



Xl + X2 



(6.5) 



where m is the lower of the constituent masses. It follows then that any multi-constituent state with a given total x 
and q± must have an invariant mass which is greater than or equal to that of the state containing just two constituents 
of lowest mass m and having the same total g"*" and x. 

We set constituent momentum cutoffs such that for states of center of mass momentum {P^ , P± = 0) the cutoff 
momenta first appear in states when the {mass)^ is 2A^ and then evaluate the range of center of mass momenta 
(P+, P±) for which cutoff momentum constituents first appear at a {raassY of A^ or higher. Thus all states with 
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(mass)^ < 2A^ and total momentum {P^ , 0) are kept; and we must then restrict the range of (P"*", -Pi) so as to ensure 
that no state with this total momentum and (jnass)'^ < A^ has a constituent beyond the cutoff boundary. This gives 
us a large domain to test whether covariance can be restored to good approximation with appropriate counterterms. 
By our theorem, if a constituent of momentum {k^ , ki±) and mass mi appears in a state of center of mass 
momentum (Pg ,0), then the state of minimum total mass containing this constituent contains a second constituent 
with momentum (^2 , fei), where /s2 = Pq — k^ and ki± = —k2±, and mass m, the minimum of the constituent 
masses. The mass of this state is 

(massY = — ^ ii H i^ (6.6) 

Xl 1 — Xl 



k 



+ 



with Xl = -j^, X2 = -^- We therefore need to characterize the boundary curve in xi, kij_ space on which 

™2 , u2 ™2 I u2 _2 ™2 i„2 

2A2. (6.7) 



ml + fcjj^ TO^ + klj_ 

Xl 1 — Xl 
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Xl 1-.T1 a;i(l-xi) 


Solving this equation for klj_ we have 




kl^ = 2A^xi{l- 


' Xl) — mi(l — Xl) — m^xi. 



(6.8) 

For A^ ^ m,l, we find klj_ achieves its maximum value very close to — , when xi = ^- To ensure no deterioration in 
the cutoff for momenta P"*" > Pq , we define the cutoff boundary curve to be 



while for xi < i 



A^ A;+ 1 

kl^ = m^ for Xl — -^ > - (including xi > 1), (6.9) 

2 Pq 2 



kj^ = 2A2xi (1 - Xl) - m^, (6.10) 



which applies only in the domain xi > ^^5- + O(^) for which klj_ > 0. We will ignore corrections to the bound on 
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X that are O(^), but a careful analysis must retain these, xi < j^ is not allowed. Note that xi is defined in terms 

of Po+, even for P+ ^ P+ . 

The reader may at first feel somewhat uneasy about the statement "including xi > 1." Consider for simplicity a 

two-body system of equal-mass particles. First consider the frame P+ = Pq^ and P± — 0. The internal mass of the 

state is 

M=^(^ + !^)P^ (Ol) 

Define xi = -^ and X2 = -j^. Thus in the frame P+ = Pq^, < xi < 1, where we have dropped terms of order 

TO^/A^. Next consider the frame P+ = P+. Let c = ■^. Thus in this frame < xi < c and c is greater than 1 

if P+ > Pq^. Wc also have X2 = c — xi and < X2 < c. If we did not extend the boundary as exhibited in (|6.S|), 
however, the x^'s would be limited as < x^ < 1; and one could not have a two constituent state with P+ > 2Pq'". 
Such a cutoff would likely be too restrictive. 

We next determine the boundary curve in P+ and P± space for which the cutoff mass is A^, namely, constituents 
on the cutoff boundary occur only in states of (mass)^ = A^ or higher. If the center of m ass momentum of a state 



is (P+, Pj_) and we consider one constituent {ki , ki±) on the cutoff boundary (6^ — 6.1C ) with xi = ki / Pq , then 



from our theorem above we know that the invariant mass of this state must be at least 

^^2 _ m^ + (fci± - £.lPl^^ , m^ + (fc2± - ^2Pi_? _ m^ + (fc^ - ^iP^)^ 

6 + 6 " ad-ci) ' ^ ^ ^ 

where ^1 = kf /P~^, £,2 = ^2 /P'^, and the second equality follows from the relations k^ = P'^ — k^ and k2i_ = P± — ki±. 
We ask, what must (P"*", P±) be to ensure that the minimum of this mass over all allowed values of ki , with klj_ on 
the boundary, is A^? 

26 



First , not e that M^ is a minimum when ki± and Pi are in the same direction. Let 7 = P~^ /P^ , so that ^1 — 0:1/7. 
Using ( 6.12 ) and the requirement that M^ be greater than A^, we find 



xjPl ^ 2xi-f \ h^ I \P±\>A^xi{-f-xi)~j^{kl^+m^). (6.13) 

Dropping the solution which excludes P± = 0, wc have 

<| F_i |< - [7 I fci± I - (A'^i(7 - ^1) - I'm')'^'] . (6.14) 



Substituting for fci^ from ( |6.9 -6.10), dropping terms of order rn^ /A'^, and minimizing the right-hand side with respect 
to xi, we find 

0<|Pj_ |< A[(l-7)(27-l)]^/^ for 1/2 < 7 < 3/4. (6.15) 

and 

0<|P±|<^ for 7>3/4, (6.16) 

2v2 

Using the definition of 7, this implies that only states with P+ > P(l /2 are allowed. 

C. Cutoffs and field operators 

To summarize the results of the previous section, the cutoff boundary on each constituent's momenta {kf, ki±) is 
specified by 

fc^=2A^^(l-^)-m^ for ^P+ < fc+ < ^ (6.17) 

and 

kf^ = ^-m' for ^<fc+, (6.18) 

where m is the lowest of the constituent masses. Then, in a considerable range of center of mass domain {P± = 
0, P+ = -|- to Pi = -A=, P+ > -^), when a cutoff constituent appears in a state, the internal mass of the state is 

guaranteed to be at least A^. This domain for physics with cutoff (mass)'^ > A^ is shown in Fig. 4. 

At all stages of the calculation we will also want to provide a buffer zone outside this constituent momentum cutoff 
boundary extending kf^ roughly a factor of 2 and kf' a factor of 1/2. We can accomplish this by setting the outside 
of the buffer zone at 

fcL=4A^^(l-^)-m^ for ^P+ < k^ < ^ (6.19) 

and 

P+ 

fc,L=A2-?7i2 for -^<k+. (6.20) 

The buffer zone allows the use of a smoothing cutoff in order to let interactions die gradually. In (6.17-6.20), the 
cutoff A eliminates both ultraviolet transverse and infrared longitudinal degrees of freedom. In order to define such a 
constituent cutoff which still limits the invariant mass of a state, however, we are forced to introduce a longitudinal 
momentum scale Pq . Dependence on this scale must also be eliminated as part of the renormalization process. 

We now have the boundary on the degrees of freedom kept in k± and fc+. This is shown in Fig. 5. These cutoff 
boundaries are to be employed in the integrals over momenta in the momentum-space expansion of the field operators. 
The quark field operator is 

^(-) = E ^ J ^KS^ Hk)e-^'- + dl{k)e^'^-] , (6.21) 
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and the gluon field operator is 

A\x) = J2 / T. , n^ k«A(fc)e-'^- + ^;:a\{k)e^^-A , (6.22) 

^ J c IbiT-^v k+ L -I 

where color indices have been suppressed. If we choose sharp momentum cutoffs, the momentum integrals may be 
explicitly written as 

dk+(fk±= / dk+ / d^k^e{2A^x{l -~x)-m^ - k]_) 

+ / dk+ d^k^0( m^-kl), (6.23) 

Jp+/2 J 2 

where now fc+„„ = ^P^. 

With sharp momentum cutoffs it is possible to develop non-analyticities in the structure of counterterms. We 
illustrate this with the following example. 

Consider the type of transverse momentum integration that occurs in the gluon self-mass. In its simplest form this 
integral is 

/= f d^ki^d^k2^S^{P^- ki^- k2^) e{A^ - kfj^) e{A^ - k^j^), (6.24) 

d^ki^ 9{A^ - kl^) e{A^ - (Pi - k2±f) , (6.25) 



2" ^' '"^^ ~ '"^ ' 2-' 
27r 



with A a very large parameter. Introducing the variable k± by ki± — k± + -P±. k2± = —k± + ^P 

1 
I = 

with (for Pi << A) 



I=^J #'«L., (6.26) 



Kmax — Minimum[A — - \ P± \ cos(j),A + - \ P± \ coscj)] . (6.27) 

Then, 

I^vergent = TtA^ - 2A \ P^ \ . (6.28) 

The sharp cutoff has given rise to a non-analytic divergence. To avoid such non-analyticities in the structure of 
counterterms, analytic cutoff functions are preferred. For example, the momentum integrals can be cutoff as 

where x — k^ /P^^ . C{y) is a cutoff function which equals 1 for y = and decreases analytically to for y = 1. The 

I 2 I 

integration over longitudinal momentum is thus cutoff at a minimum k^^^ — ^tt^q • As long as C is analytic the 
self- mass will have to be analytic in Pj^; the | P^ | divergence will disappear. 

VII. SIMILARITY RENORMALIZATION SCHEME 

The bare cutoff Hamiltonian will be solved in two stages. The first stage is a renormalization stage in which effects 
at relativistic momenta are computed using perturbation theory. These effects are computed to a specified order in g 
in perturbation theory. For the simplest computation the renormalization calculation is limited to second-order terms 
in g. To study the next level of complexity, the similarity renormalization program can be carried out to fourth order 
in g, but that is not attempted here. 

The output of the renormalization stage is an effective Hamiltonian which is dominated in weak coupling by 
nonrelativistic relative momenta. The bound states of the effective Hamiltonian in weak coupling are pure gg, pure 
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qqq, or pure gg states, exactly as predicted in the CQM. All major effects of gluon emission and absorption — or more 
complex processes — are absorbed into the effective Hamiltonian for these states. Once the effective Hamiltonian is 
constructed, the computation of bound state energies involves readily executed numerical computations. 

A new procedure for the rcnormalization of Hamiltonians has been developed by Glazek and Wilson to best meet 
the needs of light-front computations |lj,^ql • The new procedure is needed to simplify qualitative calculations leading 
to the effective Hamiltonian. Of especial importance is an ability to handle linear potentials correctly. The problem 
with a linear potential is that even if it has a small coefficient (even of order g^), the linear potential becomes of 
order 1 or even larger simply by going to large enough distances. But once a linear potential has become large enough 
it can no longer be treated perturbatively and must instead be incorporated into an unperturbed Hamiltonian, for 
either qualitative or quantitative analysis. Since a crucial question is how linear potentials do or do not develop, a 
procedure which handles them sensibly is of critical importance. 

Another requirement for the new rcnormalization procedure is that no small energy denominators should arise in 
any order of perturbation theory. That is, energy denominators of the form 1/{E — E') should never arise with \E — E'\ 
being far smaller that either E or E' . Normal perturbative treatments of field theoretic Hamiltonians lead constantly 
to such small energy denominators. But for exploring the complex phenomena of light-front rcnormalization, which 
will require crude and qualitative analyses, it is desirable to avoid encountering the possibly large effects resulting 
from such small denominators. Such large effects are very difficult to estimate using qualitative arguments. Large 
effects arise in particular when kinematic constraints force E — E' to have a single sign yet allow it to be small. A 
typical example is coUincar emission of gluons at large transverse momenta. In this case, the incoming and outgoing 
constituents all have large light-front energies due to their large momenta. Nevertheless, the final state energy is only 
slightly above the initial state energy, yet cannot ever be lower than the initial state energy. 

The third requirement we place on the new renormalization procedure is that the effective Hamiltonian should have 
a Jacobi-type cutoff severely restricting internal momenta while the center of mass momenta are restricted only by 
the original bare cutoff. This will allow us to verify the restoration of longitudinal and transverse boost invariance in 
the effective Hamiltonian itself. 

The first stage computation also identifies all cutoff dependencies in the bare Hamiltonian, and therefore can be used 
to identify all needed counterterms to eliminate the cutoff dependence in perturbation theory, including dependencies 
that violate as well as conserve boost invariance. 

A. The similarity transformation of Hamiltonians 

For now our discussion will be valid for any Hamiltonian with non-negative eigenvalues. Specialization to the 
light-front case occurs below in Section VH.B. The basic idea of the "similarity renormalization scheme" is to develop 
a sequence of infinitesimal unitary transformations that transform the initial bare Hamiltonian Hg to an effective 
Hamiltonian H^-, where a is an arbitrarily chosen energy scale. When the transformations are multiplied together, 
one arrives at a transformation So- '■ 

H, = S^HbSI (7.1) 

The basic goal for the transformation So- is that Ha^ should be in band-diagonal form relative to the scale a. What 
this means, qualitatively, is that matrix elements of iJ^- involving energy jumps much larger than a (other than jumps 
between two large but nearby energies) will all be zero, while matrix elements involving smaller jumps or two nearby 
energies remain in H^. 

We require that for cr — > oo, Ha^ — > Hb and S'o- — > 1. The effective Hamiltonian we seek involves H^ with a of order 
the quark or gluon mass. The similarity operator S^ may be expressed in terms of the anti-Hermitian generator of 
infinitesimal changes of scale To- as 

/■OO 

S^=Texp{ da'T,,), (7.2) 

J a 

where T puts operators in order of increasing a. Clearly, this satisfies the above limit condition. Using 

T„ = S,^^-^Sl (7.3) 



the differential form of (7.1) is 



^ = [H.,T,], (7.4) 
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which is subject to the boundary condition hnig-^oo H^ = Hb from above. We shall never need to explicitly construct 
the similarity operator Sa in what follows. 

We need to specify the action of T„. We write Hb = Hq + Hf, Her = Hq + Hj^, let Ei and Ej be eigenvalues of 
Hq, and denote any given matrix element of H^ as H^ij. The differential equations of the similarity renormalization 
scheme |l5[ that determine both H^j and T^j are, 



dH„ 



da 



T — 

"■'■' ~ E, 



1 



da 



(}nfcnj)H„ij, 



E,, 



(1 - faij)[Hla,Ta]ij - —{In faij)Hcr,,j 



(7.5) 



where faij = f{xai]) with 



i-o-ij 



\E.-E,\ 

E; + Ej+a' 



(7.6) 



These equations are consistent with (7^). The function f{x) should be chosen as follows: 



0<x<i, f{x) = l 



(near diagonal region); 



1 2 

— < X < —, f{x) drops from 1 to [transition region); 



3<x<l, f{x)=0 



(far off diagonal region) 



(7.7) 



f{x) is to be infinitely differentiable throughout < a^ < 1, including the transition points 1/3 and 2/3. 

It will be shown that with these definitions -ffo-ij is zero in the far off-diagonal region and T^ij is zero in the near- 
diagonal region. The second claim is immediately obvious from the form of To-y in (7.5), for 1—/ and do- (In/) vanish 
when < a: < 1/3. To see that H^ij vanishes when x > 2/3, we first note that since Hq is diagonal, d^ihi f)HQij 
vanishes identically (/ is 1 and d^f is for x = 0); and from this result we can indeed declare that the unperturbed 
part Ho of H„ does not vary with a, with T„ being of order Hi or higher. Now we can re- write the equations (7.S) as 



r„ 



rf^lA^^'"''^"'^''''^"^'" 
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E, - E, 



(1 — I<yi])[Hlc^Tc\ 



, djaij , 1 ^ 

da fcrij 



We will see below that Hj/f remains finite even when / = 0. Solving the first equation with the stated boundary 
condition at cr — > oo gives 



H, 



laij 



= /< 



(TIJ 



ji/f?^-^ da'[iI/,.,T,,],,| 



(7.9) 



Since f{x) vanishes when \x\ > 2/3, we find that Hdj does indeed vanish in the far off-diagonal region 



That T„ij is zero in the near-diagonal region means that a perturbative solution to H, 



(71] 



in terms of Hf^; will never 

involve vanishing energy denominators, which are a potential source of large errors in other perturbative renormal- 
ization schemes. That H„ij is zero in the far off-diagonal region will help identify divergent terms and determine the 
form of counterterms necessary to remove these divergences. In LFQCD, we are considering cutoffs which severely 
break the symmetries of the canonical theory. The renormalization process is therefore sure to be quite complicated, 
and the advantage of the similarity renormalization procedure is that it allows a clean identification and separation 
of divergences. 

Finally, we have to discuss renormalization conditions. The bare LFQCD Hamiltonian has a limited structure 
— canonical terms, an artificial potential, and renormalization terms. The renormalization counterterms are to 
be determined by fixing outputs: both physical masses and renormalized coupling constants defined in terms of 
measurable parameters, and a subset of the constraints of Lorentz covariance. These outputs are to be obtained 
perturbatively by solving the Hamiltonian bound state problem. But for now we mainly want the functional form of 
the counterterms, which has to precede determination of their strengths. The main requirement for counterterms in 
Hb is that they remove divergences in the integral over the commutator [Hiui ^Tc']. 
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We can summarize the equations for Hj^ and To- as 

Hi„ = hZ + [Hi^,,T„,] , (7.10) 



where Hf^^, — faijHf^, and the Unear operation R is 




da'X^'y. (7.11) 

I Jo- 

Using the equation for Hju ahows us to write 

T„=Hf^T + [Hia'.T„'] , (7.12) 



where 



and 



^^-^^^ 'K—^, [d^^^'' I ^'^' ^^-^^^ 



^0 



d . \ r , ,^. 1 



/,,, / da'X,,,, + (1 - fa^J)X,,J. (7.14) 



w:^^ Ej - E, \da-'"''J J„ " '■' E, - E 

This equation structure saves us from writing the solution for T^ separately because it is obtained by substitution 
from the equation for Hj^r; namely, Hf^ -^ ^EtT' ^^^ ^^^ higher order terms have the substitution ... . So we 

can start writing the iterated solution for Hi„: 

Hi, = HP, + [HP,, , hP,,t] (7.15) 

^ « 'r 

+ [[Hi,»,Hj,„j.] ,Hj,,j,] , 
^ 'r' 

' 'R 

+ [HP,„[Hp,,„HP,„r] ] +... 

. 't' 

' 'r 

The counterterms in HP must be chosen to cancel the divergences which occur in integrals over intermediate states 
at higher order in HP , and such counterterms must also then be included in higher-order iterations. If a limit to this 
process exists, the Hamiltonian is said to be rcnormalizable. 

B. Discussion of the scheme 

Now we may discuss the similarity rcnormalization scheme in more detail. The ultimate aim of this scheme is 
to transform the Hamiltonian Hb into a manageable, band-diagonal form H,. Hb is the bare cutoff Hamiltonian, 
forced to be finite by the imposition of some cutoff A. The differential transformation framework produces a set of 
Hamiltonians H,,, where a' ranges from ©(A) down to some scale u and thereby dresses the Hamiltonian. One step 
in this process is the determination of the form of the counterterms which must be included in Hb so that each of 
the matrix elements of the transformed Hamiltonian H,ij has no large dependence on A. This is the rcnormalization 
process: as we send A -^ oo, we get a renormalized, scale-dependent effective Hamiltonian H^. This does not finish 
the rcnormalization of the Hamiltonian, however, for the finite parts of the counterterms in Hb will produce in H^ 
unknown constants and functions of momenta which must be adjusted to reproduce physical obscrvables and to restore 
the symmetries which were broken by the cutoff A. These quantities are to be fixed, then, by solving H^\ and one 
should be able to do this with a combination of few-body Hamiltonian methods and weak-coupling diagrams. 

We show this rcnormalization scheme pictorially in Fig. 6. The grey area is the region where a is large enough so 
that H, is equivalent to the bare Hamiltonian Hb- The necessary condition is that \x,ij\ < 1/3 for all i and j, which 
is true for 
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where A is the maximum and jjl is the minimum energy of the free states of the cutoff Hamiltonian Hb- By lowering 
the scale a below this region for a given A, we start transforming Hb — eventually producing a band-diagonal 
Hamiltonian at some sufficiently low scale. Thus we form a new "Triangle of Renormalization," as shown in Fig. 
6. First, we eliminate degrees of freedom by fixing the cutoff at some value Ai, and we denote the corresponding 
bare Hamiltonian as Hg . Now we perform the unitary similarity transformation to bring the Hamiltonian to a band- 
diagonal form characterized by some scale ctq, which we denote as iJg in the figure. Next, we increase the cutoff to 
some value A2 to get H^ and then transform via a to get iJg, defined at the same similarity scale ctq as was i?Q. 
As we increase Am in this way, we end up with a sequence of band-diagonal Hamiltonians H^ . The renormalized 
Hamiltonian at scale ctq is the limit of this sequence, H^ = limjv^oo ^o^- 

Of course, we could use the similarity transformation to bring the effective Hamiltonian to any scale (T„ and obtain 
a set of sequences {H^}. As we change a, we change the characteristic scale; but the physics is invariant with respect 
to this change for large enough Ajy. The utility of the similarity renormalization scheme is that the transformation 
is invertible, as indicated by the double-arrowed lines in Fig. 6. Thus, if we find a Hamiltonian that is finite and A- 
independent for any one scale am the differential similarity framework guarantees that we can obtain a Hamiltonian 
that is finite and cutoff independent for all cr. Note that we require that each matrix element of Ha be cutoff 
independent for external momenta which are small in comparison to the cutoff, which is more restrictive than just 
requiring this of the eigenvalues of H^. 

We now discuss the choice of counterterms. We start with the bare cutoff Hamiltonian HB^ which consists of the 
canonical Hamiltonian with cutoff A plus a set of counterterms which have an explicit A-dcpendcnce arranged so that 
physical observables obtained from Hb are A-independent. These counterterms are at the outset undetermined. As 
we perform the similarity transformation, we get an effective Hamiltonian H^, where the scale a replaces A. This 
is done through a precise differential framework. Qualitatively, we are separating high-energy degrees of freedom (of 
order A) from low-energy degrees of freedom (of order a). In the language of the phase space cell analysis of Section 

V, we have A ~ transverse (longitudinal) widths {5y\^^^{5y~) respectively and a ^ transverse (longitudinal) widths 
{5x±)~^{5x~) respectively. Thus by varying A, we determine which divergences have to be eliminated from H^- The 
counterterms which eliminate these divergences must be added to Hb, and so can have no dependence on a. The 
structure of each counterterm H'j^f, is limited by the structure of the divergence it must cancel — namely, by the 
external legs associated with states i and j. Along with the diverging part of each counterterm, we may in principle 
associate an unknown finite piece, whose precise form is to be determined by fitting data and restoring symmetries. 
These terms must be restricted so that a piece which removes an ultraviolet divergence, for example, will have a precise 
behavior with respect to external transverse momenta, but may have a finite piece associated with it which contains 
an unknown function of the external longitudinal momenta — provided that function does not break boost invariance. 
A similar rule applies to counterterms for infrared divergences. However, remember that the individual constituents 
associated with external momenta do not correspond to physical states, and so we cannot include counterterms for 
apparent infrared divergences which would be eliminated upon integration with a test function. 

Finally, we need to fill in some details in order to match the general presentation of the similarity scheme given 
here to the particular choices of light-front coordinates and the cutoffs on constituent momenta described in Section 

VI. The bare LFQCD Hamiltonian has a cutoff A which has the dimensions of transverse momentum. The cutoff 
A of Section VI is designed to simultaneously eliminate constituents with large light-front energies due to either 
large transverse or small longitudinal momenta. We want to define the light-front similarity scale to have the same 
dimensions as A and to be a true invariant mass scale. The scale a is defined through (7.6), which we now re-define 
for the light-front as 

Here, P~ and P~ are the (off-mass-shell) light-front free energies of the states i and j (namely, sums over the states' 
constituents' light- front free energies); and 

_2p2 I ^2 

p- ^ 'y " . (7.18) 

Since {P^ ,P±) is conserved in any process — that is, Hij ex S{P,j^ — Pj')5'^{Pii_ — Pji_) — we see that Xaij as given 
in (7.18) is independent of the total longitudinal momentum P^ and transverse momentum P^. Defined in this way, 
the light-front similarity scale a' interpolates from a mass scale C(A) down to some final scale cr, w hich s ets the mass 



scale of the bound states of the effective Hamiltonian. Moreover, from (6.1) and (6.2) we see that ( 7.17 ) reduces to 



\M]-Mf\ 
Mf + hq + a2 ' ^^-^^^ 
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with Mi and Mj the total masses of the states i and j. So the similarity transformation replaces the cutoff A, which 
explicitly breaks boost invariance, with a boost invariant mass scale a. 



VIII. EXAMPLE CALCULATION: SECOND-ORDER GLUON MASS CORRECTION 



In this section we provide an example of the application of the new regulator scheme and the similarity transfor- 
mation perturbation theory to second order in the coupling constant. We calculate the correction to the gluon mass 
coming from intermediate two-gluon states, which is d epict ed in Fig. 7. This correction contains severe divergences. 

First, we need to find the explicit form for ij'^' in (7.15). After substitution from the previous equations and use 
of the Hermiticity of Hf , we find 



^laij — [Hla'iHicr'rVj — / , ^lik^lkj 
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9<yijk 
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Pu 



'R k l ^ « ^ k ^ j k ) 

where the sum over intermediate states k includes an integration over momenta. The function g is 

d 



9<yijk — /< 



ai] 



da fa'ik-r^fa'jk, 



(8.1) 



■2) 



and it serves to keep i and j in the near-diagonal and transition regions {x^ij < 2/3) and j and k in the transition 
and far off-diagonal regions (xg-jfe > 1/3). The identity fa-ij = faji and integration by parts give 



Qaijk I 9<7jik Jaij\^ Jaikjajkjt 
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SO that we can write 
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(8.4) 
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2iPr^p-)(P7 
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aijk 



9(7 jik)- 



The first term is the same as that obtained in the Bloch rc nor malization scheme ITS], but now modified so that the 
energy denominators never get small. The second term in ( ^.4[ ) vanishes as P[' —^ P~ . Divergences occur when P^ 



is much greater than P^ and P- , in which case the first term in (3.4) will give the dominant contribution to (8.1). 
For the gluon mass correction, the initial and final states have the same momenta, so that 



TjC^) -\^ IjB TT 



1-/' 



p: 



Ikj p- 



(8.5) 



We denote the single gluon state as \i) = |P, a. A), where P = (P+, Pi), a, and A are its momentum, color index, and 
helicity index. Using the diagrammatic rules for the canonical vertices M, we have: 
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9^Ca 



2P+16 



'—S'^" j dk+d^k^^dk + d?k2^5\P - fci - k2)M2XX' {P. ki,k2){l - f^ik) 

'' J C 



(8.6) 



where the factor ^ is the symmetry factor, P- ~ ™" p+ ^ ' -^fc =^1+^2 ^ "'^ + ^^ + "''^ + ^^ , /aifc = f{x^ik) 

with / given by (7.7) and x^ik given by (7.17-19), CaS"'' = f'^'^f'^'i = Nd"^, and M2\\' denotes the two vertices 
with the energy denominator. 

To find the mass correction, we set P± = and define 



rr(2) I 



^^AV^-, 



where the mass shift is Stuq — Sitiqj^ + Sttiq^ . 

Integration over fc2 and the cutoff scheme of Section VI lead to 
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and the cutoff function 0^{k) is defined via (6.23) as 
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(8.10) 



with e ~ ~2A^' -^ ^^^ bare cutoff, and m the smaller of mc and ttif. We use sharp cutoffs here only for the purpose 
of illustration. We stress that smooth cutoffs are preferred so as to avoid non-analyticitics. 
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(8.13) 



So we see that the external (total) momentum P+ enters only through the cutoffs. The term independent of fa-ik 
in (8.12) gives rise to divergent contributions as the cutoff' A becomes large, as well as finite contributions, and the 
second term involving /^^j, is a finite, cr-dependent contribution (fa-ik cuts off the integration before it reaches the 
bare cutoff A). 
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Evaluating / for the case P"*" > P^ , for example, the divergent part of the mass shift is 

6ml^ ^-N£^{Fi+F2+F3+Fi). 



Fi is the quadratic ultraviolet divergent part, 
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F4 is the pure infrared divergent part. 
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Note that the mass shift is negative, as it must be in second-order perturbation theory. 

The coefficients of the ultraviolet quadratic and logarithmic divergences are multiplied by functions of longitudinal 
momenta which are in accord with the power counting rules but which are not boost invariant. We have to subtract 
them entirely, and the finite parts of these counterterms are just arbitrary constants since an arbitrary function 
of longitudinal momenta violates longitudinal boost invariance. The divergences arising from small longitudinal 
momenta also violate longitudinal boost invariance because of the cutoff scheme, and so these divergences are also to 
be subtracted away, leaving room for only an arbitrary constant. 

The evaluation of the second term in (8.12) depends explicitly on the choice of the function /. For qualitative 
purposes, we may choose a step function, namely. 



which results in the effective cutoffs for the Jacobi momenta 



(8.19) 
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These yield a result bounded by a and finite as A — > 00. 

Thus, with the cutoff theory, even if we set the bare gluon mass to zero, the gluon mass correction Sitiq does not 
vanish. After subtracting the infinite mass correction, there remains a finite contribution which depends on the mass 
scale parameter a. This dependence is altered by spectators because Xa-ij changes when spectator energies change. 

IX. INFRARED COUNTERTERMS IN LOWEST-ORDER PERTURBATION THEORY 

In Section IV we discussed the infrared counterterms in the canonical Hamiltonian itself. Products of the interaction 
Hamiltonian lead to light-front infrared (small fc+) divergences in addition to the ultraviolet (large k±) divergences. 
Note that we do not have the conventional infrared divergence (small k±) since gluons are massive. As we have 
pointed out there is a crucial difference between these infrared divergences and ultraviolet divergences. 

In constructing counterterms for the infrared divergences in the products of the interaction Hamiltonian, we find 
strong cutoff dependence already in tree level amplitudes — for example, in the 0{g'^) quark-antiquark scattering 
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amplitude. Here, what one has in fact are infrared singularities. When we diagonahze the Hamiltonian, these 
amphtudes are integrated with wave packet functions, and nonintegrable singularities in the amplitudes give rise 
to energy divergences. Thus to find the infrared counterterm it is not enough just to isolate the infrared singular 
amplitudes; we have to integrate the amplitude over the external leg with a wave packet function. For QED, in 
contrast, where the external legs correspond to physical particles, all tree level singularities have to cancel. No such 
cancellations will be assumed here. 

In this section we outline the construction of infrared (small k^) counterterms up to order g^. The infrared 
divergences arise from cither an exchanged infrared gluon or an exchanged infrared fermion in the intermediate 
state. The counterterms arising from infrared gluons are of interest when looking for the origin of confinement, and 
the counterterms arising from infrared fermions are of interest when looking for the origin of spontaneous chiral 
symmetry breaking. 

A. Counterterms from infrared gluons 

1. Counterterms in the quarle-antiquarle sector 

The interaction Hamiltonian Hqqg gives rise to an effective Hamiltonian in second order whose matrix element for 
the qq states \i) — |pi, si, ai;p2, — S2, 02) and \j) = |p3, S3, a3;p4, — S4, 04) (see Fig. 8) is given by 



H\i{^= -g' r„"3„,r„-,„,M2., ^2 U.J (1 - U^kU,k) + o((P,7)-3) 



(9.1) 



where pnjSnj oin ^i"C the nioniciituni, helicity and color index of the nth quark, 
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The function 9^(k), which arises from the cutoffs on the intermediate gluon, is defined in Section VII: 
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r-, A the bare cutoff, and m the smaller of mc and mp. Here we have also defined 
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Pj, denotes the intermediate state energy: 



Pf^ = p^ + q + P2 for pt > pt and P^ = p^ - q + p^ for p+ < p^ 



(9.7) 



36 



The expression for Ti can be further simphfied, 
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LPi -P3 -9 P4 --P2 

To identify the infrared counterterms, we look for singular terms. We have, explicitly, 

K ^Vx^ V2 , Pi = Ps + Pi ■ (9-9) 

The intermediate state energy is dominated by that of the infrared gluon, namely, 

FjT « g" for p+ > p+ , PjT « -q" forp+ >p+. (9.10) 

As the gluon longitudinal momentum fraction reaches e, both Xaik and x^jk become greater than g; and as a 
consequence 1 — faikfajk -^ 1- The leading singular term is 

2 f 

.2j,a J^a g g (^ ^]]:G A Mj 

4ff Ja3ai^04a.''«i^3«s.s.^i (p^^ _ P3^)2 + ^2^ j (p+ _ p+)2 

X 9^{p,-P3)+0Hp3-Pi)]- (9.11) 

To find the infrared divergence we integrate over p3 with a test function 4>{p^ , P3±) , which we may assume is nonzero 

L. This evaluation is complicated by the cul 
(9.5), we find a divergent piece 
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only in the range where jaij — 1. This evaluation is complicated by the cutoff function. Using the definition of G\ 
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where k — p^ —pi; and a similar piece for p^ —p^ > e, which is the same as (9.12) except that (t){p + k) -^ (j){p— k). To 
find the small-e dependence in I, we write fc+ = ex, switch the order of integration, and expand in e. The divergence 
has the form 



2m2 /■ .T. fc2 



'^'^^7I^^::Wl^^ITT<^(P^PiJ- + ^^)- (9-13) 



Hence we need to include a counterterm for this linear divergence in Hb'- 

^9'T^,.J:..Js.sJs.s. — -. (Pi^-P^^)' M-Pt). (9.14) 

Without our cutoffs and for zero gluon mass we would find exactly the same counterterm but with the opposite sign 
as the counterterm for the instantaneous four-fermion interaction in the canonical Hamiltonian, and hence no 0{g^) 
counterterm would be necessary for the linear infrared divergence. But our cutoff scheme, which is dependent on the 
choice of massive constituents and thereby eliminates zero modes, does not allow complete cancellation to occur. 

There are also -^ type singularities which lead to logarithmic infrared divergences. But these divergences are 
cancelled from the two ^-functions in (9.8). Thus there is no counterterm necessary to remove a logarithmic infrared 
divergence in qq scattering to second order |Q. However, as discussed in Section IV. B for the canonical Hamiltonian 
counterterms, the use of a symmetric cutoff to ensure the cancellation of divergences from small positive and neg- 
ative longitudinal momenta excludes any possibility of finite contributions from exact zero modes. To counter this 
elimination of zero modes, we can include a finite term analogous to (4.17). Here the transverse structure in M2 — 
for example, the product of the first and last terms in (9.2) — may keep such terms from vanishing, in contrast to 
terms corresponding to instantaneous zero-mode gluon exchange. The finite counterterms also involve a product of 
two fermion color charges as one needs to start building true confining potentials. True confining potentials confine 
only non-zero color charge states as opposed to the artificial potential which confines all states. 
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2. Counterterms in the quark-gluon sector 
Now we consider the qg states. The effective Hamiltonian in second order (see Fig. 9) is given by 
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Again, to identify the infrared counterterms, we look for singular terms. As in the case of the quark antiquark 
sector, the intermediate state energy is dominated by the infrared gluon and the factor (1 — faikfajk) is replaced by 
1. The most singular term (leaving out the factor fa-ij) is 
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To find the divergence we must again integrate over k2 with a test function, which again yields a linearly divergent 
term as above. Here also the cutoffs and nonzero gluon mass spoil any complete cancellation of this divergence with 
the counterterm for the two quark-two gluon instantaneous interaction. 

After the subtraction of the linear divergence, a linearly singular term (proportional to vtiq) survives; but as in 
the case of the qq scattering amplitude, there is no divergence after integrating the amplitude with a test function. 
There are also other l/g"*" terms in (9.18), but the logarithmic divergences they produce are cancelled from the two 
0^-functions. However, as with the qq sector, even though the logarithmic divergence for small q^ is cancelled by 
that for small — g^, we need to include a counterterm to account for the possibility of finite effects from exchange of 
gluons with exactly (j+ = 0. 



B. Counterterms from infrared fermions 



In this subsection, we consider the contributions to the matrix elements of the effective Hamiltonian in the quark- 
gluon sector arising from quark exchange. First, we discuss the contribution from an intermediate one-quark-two-gluon 
state, as is shown in Fig. 10a. The matrix element is given by 
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with P3 =Pi - kj, p3±_ = pi±_ 
In the limit pg ^ 0, we have 
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which leads to a logarithmic divergence when we integrate the effective Hamiltonian with wave packet functions which 
are functions of external momenta. 

Next we consider the contributions to the matrix elements of the effective Hamiltonian in the quark-gluon sector 
arising from antiquark exchange, that is, from an intermediate two-quark-one-antiquark state, as is shown in Fig. 
10b. The matrix element is given by 
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with p+ = fc^ - p+, p3^ = k2i_ - pii_ 
In the limit pg -^ 0, we have 
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which leads to a logarithmic divergence when we integrate the effective Hamiltonian with wave packet functions which 
are functions of external momenta. 

Note that if we keep both types of intermediate states, the singularity acquires a principal value prescription since 
the ^-functions sum to one, and thus no counterterm is needed. This is because the singularity occurs when the 
intermediate state energy is dominated by the exchanged quark or antiquark, in which case it does not matter if the 
other particles in the intermediate state are quarks or gluons. However, it is clear that the finite part of these diagrams 
will be affected differently, for the function 1 — faikfajk in (9.23) will now cut off the two diagrams differently since 
the intermediate state energies are different when the exchanged quark and antiquark have finite energy. Thus we 
expect that at higher order there will be diagrams which contain quark or antiquark exchange of finite energy and 
have divergences associated with another process. Counterterms for these diagrams will have a different dependence 
for the intermediate qgg and qqq states and so a logarithmic divergence will not acquire a principal value prescription 
as here in second order. The counterterms necessary to remove such higher-order divergences will provide a possible 
source for explicit chiral symmetry breaking terms in the renormalized Hamiltonian. 

Moreover, as with gluon exchange above, despite the cancellation of near-zero k^ logarithmic divergences, we must 
include a finite term which might restore any physics lost by the elimination of exact zero modes. As discussed 
in Section IV. B, there are finite non-canonical terms which have the structure of these logarithmic divergences and 
which satisfy the requirements of power counting and boost invariance yet break chiral symmetry. Since such terms 
necessarily contain the effects normally associated with the vacuum, it is no surprise that they provide a source for 
chiral symmetry breaking (see the discussion of the sigma model in Appendix A) . 

Finally, a similar result applies for the matrix element describing quark-antiquark annihilation into two gluons, 
which involves the quark exchange shown in Fig. 11a and a similar diagram with an intermediate antiquark. We 
should here include a finite term so as to counter the removal of zero-mode fermion exchange, and this may include a 
piece which explicitly breaks chiral symmetry. It follows that the chiral symmetry violating term in the qqg coupling 
will have an arbitrary size due to renormalization from such a qqgg counterterm, as is depicted in Fig. lib. This 
term then need not vanish in the limit of only spontaneous chiral symmetry breaking. 
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C. Many-body infrared counterterms 

While we provide no thorough discussion of higher-order counterterms, we need to clarify one qualitatively new 
feature that arises beyond second order. In Section IX. A we showed that the matrix elements for gluon exchange 
include pieces that diverge like l/(q^)^, where q~^ is the longitudinal momentum exchange. This leads to divergences 
even at tree level. These divergences require counterterms, and we have argued that the finite part of such counterterms 
may contain confining interactions unless boost invariance precludes this. 

One example should adequately illustrate what happens at higher orders. In Fig. 12a we show a fourth order tree 
diagram in which two gluons are exchanged between three quarks. To evaluate the matrix elements that arise in the 
effective Hamiltonian we must integrate over small q"*". Keeping only the most singular parts of the matrix elements, 
we encounter longitudinal integrals of the form 



qt Je qi \q{ J \q2 J ' ~ \qi qi 

The entire product of matrix elements is relatively complicated, but this integral is easily analyzed. There is a 
log^leMl/M"^) divergence, and a log^ {eM2 / M^) divergence. There is actually a nested divergence in this diagram, 
which would be subtracted if one adds the appropriate third order diagram; however, even after this subtraction one 
is left with log^{e) divergences. 

Another simple example is shown in Fig. 12b, in which there is a gluon exchanged between two quarks while one 
quark interacts through an instantaneous potential with a spectator. Regardless of the form of the instantaneous 
interaction, this diagram diverges logarithmically unless a symmetric cutoff on longitudinal momenta is used. As 
argued above, the presence of this potential logarithmic divergence indicates that the elimination of zero modes may 
produce an interaction with the operator structure in Figure lib. This three-body interaction may be long range in 
the transverse direction, and it may be as strong as the two-body confining interactions. As we pointed out earlier, 
such long-range many-body interactions are required to cancel unphysical van der Waals forces if there is a confining 
two-body interaction Q . 

X. LFQCD BOUND STATE COMPUTATIONS 

Consider how bound states are handled in QED p6[ . There one includes in the unperturbed Hamiltonian an 
instantaneous potential whose form is obtained from the general two-fermion four-point function. The precise choice 
of the potential to include in the unperturbed part is a matter of some art, since one must weigh the competing 
demands of making the unperturbed calculation as simple as possible while at the same time including a large part 
of the relevant physics and ensuring that perturbative corrections are small. Invariably, the bound state calculation 
reduces to the usual nonrelativistic Schrodinger equation in the appropriate limit, and the corresponding nonrelativistic 
states are used in the perturbative expansion. The precise details of any particular calculation will vary, and the most 
efficient methods can only be arrived at after some experience. 

We want to take the same approach to bound states here, but we know from phenomenological constituent quark 
models that the confining potentials we need for the bound state calculation are far different from the Coulomb-like 
potentials that arise at low order in the usual relativistic perturbation theory. We expect counterterms for the light- 
front infrared divergences to be the ultimate source of the needed potentials. For artificially small coupling, however, 
the artificial potential is the dominant confinement mechanism. The similarity renormalization scheme and massive 
constituents with cutoffs allow us to treat this confining potential in a weak coupling framework. 

In this section, we outline the construction of the effective LFQCD Hamiltonian for low-energy hadrons and discuss 
the increasing levels of complexity in the LFQCD bound state computations in this formulation. 

A. The effective LFQCD Hamiltonian and bound states 

In the similarity renormalization scheme, the effective LFQCD Hamiltonian that will be used to compute the 
hadronic bound states is 



J^aij — ^IITT-A^oo J( 



aij 



Hf.+T.lHR,HF4T^ + T^) + ■■■]■ (^«-i) 



k i k 
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H„ causes transitions only between states staying "close to the diagonal" due to the factor /o-ij, and the effects 
of transitions to and from intermediate states which are "far off-diagonal" have explicitly appeared in the effective 
Hamiltonian H„ as a perturbative expansion. 

Recall that H^; has the cutoff A which violates both longitudinal and transverse boost invariance. First we identify 
the counterterms that must be included in H^ so that the matrix elements of the effective Hamiltonian H^ij have no 
divergent dependence on A. Then we send A ^ oo so that cutoff dependence is removed from the effective Hamiltonian 
Hij. As in standard perturbative renormalization theory, A dependence will be removed order by order in g„, where 
g„ is the running coupling constant at the similarity scale in Hamiltonian matrix elements. This avoids having to 
solve the nonperturbativc bound state problem to identify and remove A dependence. 

There is a question whether H„ij satisfies boost invariance after A — > oo. There can be finite terms in H^j which 
violate boost invariance yet cannot justifiably be subtracted. Such terms might result, for example, as a byproduct 
of divergent terms of the form In -^ f(P^), where P± and P+ are external momenta. The In A^ divergence must 

be subtracted. But the finite term In P]_ f{P^) cannot be subtracted because no arguments exist that would justify 
such a non-analytic Pj_ dependence resulting from the effects of states above the cutoff. In this case violations of 
boost invariance can only disappear at special values of g where the coefficient of the boost violating terms vanish. 
Clearly one such special value has tohe g = gg, as part of the restoration of covariance at 17,5. 

In the following it is assumed that the effective Hamiltonian H^ has boost invariance simply as a result of taking 
the A ^ 00 limit. Ha is generated as the first step in solving the bare cutoff Hamiltonian Hb- The second step is to 
construct bound states from H„. We may write the bound state equation as 



y^ Hcri]1pNaj = ^NlpNcri, (10-2) 



where the sum over j includes a sum over Fock space sectors and integrals over momenta. The light-front energy of 
the Nth eigenstate is 

£m = [M^N + PniMP^, (10.3) 

where P^ and Pmi_ are the total longitudinal and transverse momenta of the eigenstate and M^r is its mass. 

We will solve this field theoretic bound state problem in the standard fashion, using bound state perturbation 
theory. This requires us to identify a part of iJ^-, -ffo-Oi which is treated nonpcrturbatively to produce bound states. 
The essential simplification that makes further calculation possible is that i/^o does not contain any interactions 
that change particle number, so that the methods of few-body quantum mechanics can be used to solve this initial 
nonperturbativc problem. All field theoretic corrections that arise from particle creation and annihilation will then be 
treated perturbatively, as in QED. The potentials that appear in H^q must be chosen to make it possible for bound 
state perturbation theory to converge; and as stated above, this requires some art even in QED. 

The bare cutoff Hamiltonian is divided as Hb = Hq + Hf . For the determination of the effective Hamiltonian H^ , 
the unperturbed part of the bare Hamiltonian H^ is chosen to be that of free massive quarks and gluons with the 
standard relativistic dispersion relation in light-front kinematics. The interaction part Hf then contains the canonical 
interaction terms Hint, renormalization counterterms Hg^, and the artificial potential Va- The artificial potential 
is itself separated as V^ = Vb -I- V^'^ , with Coulomb and linear parts, Vq = Vc + Vl, and counterterms V^^ which 
remove unphysical effects resulting from the choice of a massive gluon. For the bound state calculation, we want to 
choose the starting Hamiltonian to contain as much of the physics in as simple a form as possible. So we define 

HaOi] = .fa-ij{Hoij + ^Oi]), (10-4) 

and write the unperturbed bound state solutions as Wj.^lyj- These bound states will be pure qq, pure qqq, or pure gg, 
as in the CQM. 

The starting Hamiltonian is determined from the the Fourier transforms of 



^0 



Vo^ Jdx-d^x^dy-d^y^[ - l-j+-(x)Vc{x,y)j+''{y) 

+ j + {x)VL{x,y)j + {y)}, (10.6) 

where Vc and Vl are given by Eqs.(2.6) and (2.7), j"*"" = jq°'+jg°' is the color vector charge density, and j"*" = j^+j^ 
is the color singlet charge density. This starting Hamiltonian determines the spectrum to order g^ and produces the 
zeroth order bound states that are used to analyze radiative corrections starting at order g^. 
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The interaction part of the effective Haniiltonian can now be analyzed perturbatively, and it includes effective 
potentials that contain the effects of interactions involving far off-diagonal states as well as renormalization countert- 
erms, which are determined order by order. Consider the effective Hamiltonian Hu generated to second order in the 
renormalizcd coupling constant. The counterterms H'jp" up to this order have been discussed in the Sections IV. B, 
VIII, and IX. They are the counterterms for the canonical instantaneous interactions, which can be combined as 



FgJ = -^ Id^ 



'B2 
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x^(^f dx-j+^{x-,x^)y, (10.7) 



and the counterterms from one-gluon exchange given in momentum space in (9.14) and (9.19), plus a similar term in 
the gluon-gluon sector. There are also quark and gluon mass counterterms, some of which we have given explicitly in 
Section VIII. The effect of these counterterms is to completely cancel the leading radiative corrections from instanta- 
neous gluon exchange, all of which diverge or produce pathologies. In addition, as has been discussed in Sections IV 
and IX, we need to include terms which counter the elimination of zero modes, like (4.24): 



Hc3 - 



g^ J dx-d^XA^dy-d^y^^\a^ ■ A^)Of{x, y){ai_ ■ A^)^. (10.8) 



Such a counterterm can break chiral symmetry, for example, if Op ^ inp (ci ■ d±). At this order. Op must be 
determined phenomenologically; it remains to be seen whether coupling coherence and/or the restoration of Lorentz 
covariance will completely determine it at higher orders and as one lets 5 — > ^s- As discussed in Section IX, there 
may be many more such counterterms resulting from the elimination of zero modes. We have not written them down 
explicitly, but this must be done to compute order g^ shifts in the hadron spectrum. 

What is the dependence of this effective Hamiltonian iJg- on the scale a? Suppose we restrict the states to only a qq 
pair. Then only number conserving (potential-like) terms in Ha^ contribute. Now if we let the scale a approach a very 
large number, only the original interactions in Hb survive since the factor {I — fa-ikfa-jk) -^ 0. Thus the canonical four- 
fermion interaction together with all counterterms survive, whereas the transversely smeared four-fermion interaction 
diminishes in strength. In the limit cr — > 00 we are recovering Hb which is no surprise. In this limit the effects of 
higher Fock sectors (for example qqg) can be recovered only by including them explicitly, and so it clearly becomes a 
poor approximation to include only the qq sector when a becomes too large. 

By lowering the similarity scale a, we reduce the allowed range of gluon momenta which can contribute, for example, 
to the binding of a quark and antiquark in a meson. These effects must appear elsewhere; and we see that through 
the similarity transformation they are added directly to the Hamiltonian via the second and higher order terms in 
Hb in (10.1). Thus by lowering ct, we put the bare gluon exchange effects of Hb into a qq potential in the effective 
Hamiltonian H^ perturbatively. This clearly changes the character of the bound state calculation. It changes from 
a field theoretic computation with arbitrary numbers of constituents to a computation dominated by an effective qq 
potential. If we choose the similarity scale a to be just above the hadronic mass scale, the major effects come from 
the qq sector. The resulting nonrelativistic calculation will not see the scale a in first approximation, since only states 
close to the diagonal — for which fcrij = 1 in (10.4) — will contribute. 

Consider Coulombic bound states with a Hamiltonian of the form H — ^ ^. p ^ - and the energy of the bound 

state scales like g'^. Since r scales like ^, if we add to the Hamiltonian a linear potential of the form ar and insist 
that the energy still scales like g^ wc infer that a ^ g^. 

The bound state equation for the unperturbed effective Hamiltonian H^q can be written as 



Lp+ pl+ J 47r 



^ ''mp mp '' "^ 



p+j. p'+^i 



+ (s)''[^-l5^1}* = i^*- (i°-») 



where P"*" is the center of mass longitudinal momentum. Using p+ ^ h^^ ^^'^ P'^ ^ 5^^' ^^'^ ^^^^^Y nonrelativistic 
bound state energy is 

Now as expected, the uncertainty principle shows that £ is minimized when r is of order -^, and p is of order g^ and 
hence £ is of order g'^. Numerical computations can easily give precise results. 
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B. Potential instability problem 

As we have argued in Section II. E, we add an artificial potential to our effective Haniiltonian to ensure that the 
bound state structure is similar to that of the CQM. An additional benefit of including an artificial potential is that 
it allows the removal of a potential resulting from the choice of a massive gluon which might otherwise cause an 
instability in the bound state calculations. 

We have shown in the previous sections that the introduction of a nonzero gluon mass leads to the incomplete 
cancellation of the linear interaction in the longitudinal direction between the canonical instantaneous interaction 
and the interaction arising from one-gluon exchange. The resulting interaction between the color charge densities is 
proportional to 

In the nonrelativistic limit, the transverse momentum difference in the denominator can be neglected in comparison 
to rriQ, and this part of the effective Hamiltonian reduces to a pure linear interaction in the longitudinal direction. 
This linear interaction may cause an instability in the bound states when they are expanded to include higher Fock 
space sectors. Taking the expectation value of this term in a color singlet state consisting of a quark and an antiquark, 
namely, |s2) = b^^djlO), we obtain the linear potential contribution to the energy. But the expectation value of the 
same interaction in a color singlet state consisting of a quark, an antiquark and gluon, namely, \s^) — [T"l .6*t(i^a°^|0), 

is the same as the previous case except that it is multiplied by the factor —^Ca + Cp — —h for SU{3). This results 
in an linear potential which gives a negative contribution to the energy. This negative contribution leads to a possible 
instability for many-body bound states. 

To remove the above instabilities in hadronic bound states, we need to subtract this linear interaction, as discussed 
in Section II. E. We add this subtraction to the artificial potential with the form 

^(l - ^-) j d?x^dx-dy- ]+''{x- ,x^)\x- - y-\3+'^{y-,x^) (10.12) 

Here we multiply by a coefficient (1 — g^/g^) rather than (1 — g"^ / g1) to ensure that the ordcr-g® linear artificial 
potential in (10.7) is always dominant for weak coupling. The introduction of an artificial linear potential in both 
the longitudinal and transverse directions then stabilizes the bound states for small running coupling constant g. It 
is not known yet whether this instability recurs for g near g,, and if so how to counter it. 

C. Sources of Complexity in the QCD computations 

A basic goal of this paper is to structure a sequence of computations in LFQCD with growing levels of complexity. 
The major sources of the complexity will be reviewed here. 

All computations envisioned here have two stages. The first stage is a computation of an effective Hamiltonian 
H„. We have defined cr to be a mass cutoff parameter, a mass above which states diminish in importance for further 
computations. The goal of the first stage is to compute H^ with a being just above the hadronic mass scale. In the 
simplest computation, H„ consists of a truncated version of the canonical Hamiltonian (namely Hf^) combined with 
further corrections (from the commutator in (7.10)) computed only to order g^ , as in (10.1). From this, then, bound 
state energies are to be computed only to order g'^ . This is a straightforward computation with the main concern to 
fix free parameters through a comparison with physical spectra after setting g'^ equal to g'^. The fit at this stage is 
extremely crude. 

At this point the results are analogous to those for positronium with only a Coulomb interaction. There are 
additional terms in He, coming from one gluon exchange that are order g^, but because of their dependence on r 
they do not correct the energy until order g® or higher. They produce spin and orbital splittings, but before these 
splittings can be fully computed one must determine corrections to the kinetic energy of order g^ , corrections to the 
Coulomb interaction of order 5'*, and one must allow corrections to the artificial linear potential of order g^ in Ha- 
Just as in QED, operators must be classified both on the basis of explicit powers of g multiplying them and on their 
implicit dependence on g arising from the fact that r scales like l/g^ in the bound states. At each stage one must 
consistently retain all terms in H^ that produce corrections of the same order in g after both explicit and implicit 
dependence on g is determined. This is standard bound state perturbation theory, with a linear potential added 
and forced to contribute at the same leading order as a Coulomb interaction, and with the underlying field theory 
producing additional interactions. 
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Actually there is an important subtlety here. One gluon exchange falls exponentially fast at large r because of the 
gluon mass. To obtain an order g^ spin splitting the gluon mass in the exchange interaction must be of order g^ or 
zero, not of order one. We allow the artificial potential to contain spin-dependent interactions at weak coupling that 
vanish as g — > gs, with a range governed by a mass that is order g^ . At no point do we allow the introduction of 
unphysically light constituents to obtain mass splittings, because there is no mechanism to prevent copious production 
of such light constituents in all hadrons as g ^ gg, and there is no evidence for the proliferation of light hadrons such 
constituents would necessitate even at weak coupling. 

The calculation of all the terms in Ha- required to compute binding energies to order g^ is complicated because 
of the large number of perturbative diagrams involved and the complexity of the cutoffs. However, there is a very 
interesting question that arises. The question is whether longitudinal boost invariance will be irretrievably violated 
by the fourth-order results. If the fourth-order results include logarithmic divergences due to the infrared cutoff, 
divergences which also have non-trivial dependencies on transverse momentum, then there will be a boost invariance 
violation. The reason is that the logarithm will involve a ratio of an external longitudinal momentum divided by 
the cutoff momentum. While the cutoff dependence can legitimately be subtracted, the finite part of the logarithm 
cannot, and so results in a violation of longitudinal boost invariance. This violation would be similar to the violation 
of scale invariance of canonical equal-time QCD at zero mass. The scale invariance violation is never eliminated, it 
instead results in asymptotic freedom. A violation of boost invariance would likewise be irretrievable in weak coupling, 
and could only go away at a special value of g. Clearly, this special value would have to be gs- 

In addition, if the fourth-order computations turn up such logarithmic longitudinal divergences, then there would 
have to be finite subtractions accompanying the subtractions of the logarithmic cutoff-dependent divergences. These 
finite subtractions would involve arbitrary functions of transverse momentum — functions that could conceivably be 
the source of transverse linear potentials. Certainly, these finite subtractions would be direct reflections of vacuum 
effects, which means linear potentials would not be a surprise. 

However, the most challenging problem to resolve is the computation of the effective Hamiltonian H^ beyond 
fourth-order. The major concern is the computation of higher-order corrections caused by removal of "wee parton" 
states. By "wee parton" states we mean states of high mass due to the presence of constituents in the infrared region 
and without large transverse momentum. The simplest example is a state with two equal mass constituents but one 
having a far smaller longitudinal momentum (smaller x in Feynman's language) than the other. Their transverse 
momentum will be set to p± and — p_l, respectively. If one constituent has longitudinal momentum p+, the other has 

2, 2 

momentum xp"*", with x small, then the states mass squared is about — — ^, where m is the constituent mass. The 
problem is that when pj_ is small, masses count, which means there can be major corrections associated with actual 
values for bound state masses. This is a problem once the binding energies are themselves of order the constituent 
masses, namely, when g is near its relativistic value gg. 

This challenge will not be resolved easily. It is unlikely that it can even be seriously addressed until after the 
fourth-order computations are complete, and one can start estimating higher-order terms and determine the extent 
to which strong binding energies result in large higher-order terms in H^. This challenge is further enhanced by the 
problem of avoiding instabilities in H^ itself for g near g^. 

XI. CONCLUSION 

We have presented a framework for the computation of bound states in QCD that is based on Hamiltonian methods. 
The essential ingredients of our formulation are the use of light-front coordinates, nonzero quark and gluon masses, 
severe cutoffs on constituent momenta, and the similarity renormalization scheme. We have argued that this formu- 
lation of the theory removes all barriers to a treatment of hadrons in QCD which is analogous to that of bound states 
in QED. In particular, we have shown that in our formulation it is natural to choose a starting Hamiltonian which 
contains many characteristics of the Constituent Quark Model, and therefore we expect that we may fashion this 
unperturbed Hamiltonian to model the basic physics well already at the outset. We take the view that the general 
principles of the CQM are true and propose to proceed from there. To achieve our perturbative starting point we 
must introduce artificial stabilizing and confining potentials to reproduce the necessary physics at small values of the 
coupling g, and we must extrapolate to the relativistic value of the renormalized coupling gs- At gs, the artificial 
potential must vanish, and the true dynamics will appear. The essential point to realize is that our formulation allows 
the possibility that this true dynamics may be discovered perturbatively. 

Thus the price we pay for our new formulation — renormalization problems from severe light-front infrared di- 
vergences and the breaking of both gauge and Lorentz covariance — may eventually turn out to be relatively small 
in comparison to what we gain: both confinement and chiral symmetry breaking may be studied with bound state 
perturbation theory. Moreover, the eventual restoration of Lorentz covariance — which requires the calculation of 
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physical observables so as to complete the renormalization process — will be made easier by the expectation that 
bound states will be well-approximated by nonrelativistic, few-body wave functions. The most daring aspect of this 
approach is the assignment of zero-mode contributions to counterterms for infrared divergences. There is no precedent 
for our analysis, and its worth can only be proven by obtaining satisfactory and relativistic results as g —^ gs- 

The possibility of combining phenomenology and field theory at weak coupling opens a wide range of interesting 
computations. One feature essential for the feasibility of the entire enterprise is our contention that a massive gluon 
prevents the unlimited growth of the running coupling constant in the infrared domain. For this to happen all pure 
light-front infrared divergences (due to fc^ getting small) need to be cancelled in the coupling constant renormalization 
calculation. This needs to be demonstrated. Then there is the computation of bound states and hadronic observables. 
The use of an artificial potential and massive quarks and gluons makes this possible to leading order in a perturbative 
expansion, with more and more of the true physics being included as one moves to higher orders. As discussed in the 
context of Lattice Gauge Theory by Lepage and Mackenzie |^^ , we must extrapolate to the relativistic value of the 
renormalized coupling in order to make the perturbative expansion legitimate. One can also readily conceive other 
instructive calculations which would not involve all the complications of low-energy QCD, such as the application of 
the framework developed here to the study of scalar field theory, QED, and especially heavy quark systems. 

The possibility of constructing low-energy bound states in QCD through a combination of relativistic perturbation 
theory and many-body quantum mechanics defines a start-up phase with a wide range of new and tractable calcula- 
tions. We have tried to emphasize those computations which are most immediately relevant and to define how these 
calculations can be carried out. The final demonstration of the feasibility of our approach awaits these computations. 
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APPENDIX A: SIGMA MODEL WITH AND WITHOUT ZERO MODES 

We will work only at the canonical or tree level here. The Lagrangian for the sigma model is 

C = Ud.ar + kd.Trf Via' + tt^), (A1) 



where 



Via'+7:')^-l^,'{a'+n') + ^{a'+7:'r. (A2) 



This Lagrangian has a continuous symmetry corresponding to the transformation 



a ^ a = (J cos a + n sma, 

t: —>■ TT = —a sin a + t: cos a. (-^-3) 



The extremum of the potential V is determined by 

^-ahM2 + A(a2 + 7r2)] = 0, 

OCT 

^^Tr[-n' + X{a'+TT')]^0. (A4) 

ovr 

For fi' > 0, the minimum of the potential is at cr^ + 7r^ = <^v = ^- There are an infinite number of degenerate vacuua. 
We can pick any one of them to be the true vacuum. For example, we can pick 

< I cr I >= cr„ 7^ 0, < I vr I >= 0. (A5) 
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Then the symmetry is said to be spontaneously broken. 

The conserved current associated with the symmetry of the Lagrangian is 



The associated charge is 



J^ = cra^TT - TTcf'^cr. 



Q = i I dx-(fxi_[ad+TT - nd+a]. 



From the canonical commutation relations for a and n fields, namely, 



(A6) 



(A7) 



(A8) 



we have 



[7r(a:),7r(y)] |^+^j^+=--e(x -y )S'^{x^-y^), 



[g,7r(0)]--*a(0), 
[g,a(0)]=z^(0). 



Thus we have 



<0| [Q,7r(0)] \0>=-ia,^0. 



Starting from 9^ J^ = 0, we have, 



dQ 

dx+ 



0. 



(A9) 



(AlO) 
(All) 

(A12) 



From (Alt) it is also easy to verify Goldstone's theorem and assert the existence of a massless pion. 
In the theory we have considered, zero modes are present; and as a result 

• there is spontaneous symmetry breaking, 

• the vacuum is nontrivial, < | ct(0) | >y^ 0, 

• a massless pion is predicted. 

Next we discuss the a model with the zero modes removed. To construct the effective Hamiltonian we cannot use 



the symmetry (A3) as a guide but can rely on power counting and locality. The Hamiltonian is 

P- = I dx- d^xx_[]-{d^(l) ■ d^cj) + di_TT ■ di_TT) + V]. 



(A13) 



Since the symmetry is broken only in the a sector, we insist that V is even in the Tr-field so that the symmetry under 
TT -^ — TT remains. Thus, since no inverse powers of mass are allowed, we can write 



y = \ml^'' + ^™'^' + Al0* + A2^^ + A3(/)V2 + A40' + \5Tt\ 



(A14) 



Since there are no zero modes, a term of the form Ag^ vanishes upon integration in (A13) and so is not included. 
The equations of motion are 



5+5~7r + d'i-K^ mln + 4:X2Tt^ + 2A3</«27r + 2A57r</), 
^d+d-(l) + dl(t)^ ml<f) + 4Ai03 + 2A37r20 + iX^cfP- + AgTr^. 



(A15) 



Next we construct the covariant current operator J^ . First we determine the canonical dimensions of the components 
J^,J^ and J^. Since the charge Q (= ^ J dx^d^x±J^) is dimensionless, the canonical dimension of J+ is J+ ^ 
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^ /^^\-2 ■ Since each term in 9'^J^ {— ^{d J~^ + d~^J ) — d-^ ■ J^) should have the same dimensions, the canonical 
dimensions of the other components are 

The components of J^ are to be constructed from the operators TT,(j),d'^,d~ , and d^ and constants. The constants 
are allowed to have dimensions of negative power of x± (as masses do) but no power of x^ . Since the operator ^q- 
does not appear in the canonical scalar field theory, we do not allow the operator ^^ to appear in the canonical 
current operator. We also implement the symmetry that J^ has to change sign under tt — > — tt. Then the allowed 
structure of components of J^ is 

J+= ai^+TT + a2(a+7r)0 + 03(^+0)71, 

J±^ bid±Tr + b2{d±(p)TT + 63(9^ 7r)</), 

J-= cid-n + C2(a"7r)0 + €3(8-^)71. (A17) 

Next we compute 9'' J^ and set the coefficients of terms involving derivatives which cannot be replaced by the 
equation of motion to be zero. Then we get 

02 = -C3, 03 = -C2, ^2 = -h- (A18) 

Further, from the structure of the equations of motion, we also have 

h = 2 — ' ^2 = 2 — ■ ^ ' 

Since the zero modes are dropped, we have the charge 

g = - / dx^cfx±[a2<l)d+Tr + asird+cp]. (A20) 



Insisting that Q generates the correct transformation laws ( AlC ) , we arrive at 

03 = -02. (A21) 

Without loss of generality we set a2 = 1 , so 

J+= aid+TT + <j)d+Tr - ird+cl), (A22) 

J-=cid-n + (l3d-TT-TTd-4>, (A23) 

J^= -(fli + ci)d±Tr + (j)dx_TT - 7r(9j_0. (A24) 

Further, if J+ and J^ arc to transform as components of a covariant four- vector, we also have to set ci — ai. Then 

JA" ^ (|)^^'^T - 7ra^0 + ai9'^7r. (A25) 

Since our effective Hamiltonian (even with explicit symmetry breaking) is supposed to give rise to the same physics 
as th at of spontaneous symmetry breaking, we demand that the current is conserved. Using the equations of motion 

^'^ Jp= -((/> + ai){mlT: + 4^X211^ + 2A3027r + 2A57r0) 

+7r(m|0 + 4Ai0^ + 2\j,tt'^(J) + 3A402 + AgTr^) == 0. (A26) 

Setting the coefficients of the different combinations of field operators separately to zero, we have 

mlai =0, ml-ml^ 2\zai = 0, 4A2 - 2A3 = 0, 4Ai - 2A3 == 0, 

2A5 - 3A4 + 2aiA3 = 0, -As + 4A2ai = 0. (A27) 

The solutions are 
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Ai = A2 = -As = -, A5 = Afli, A4==Aai, m'^^ml + 2Xal, (A28) 

and 

ai = or m^ = 0. (A29) 

If we take ai = 0, then the potential is reduced to the canonical form (A2), and the current is also of the canonical 
form. This corresponds to the full canonical theory with a symmetry preserving vacuum and a doublet in the spectrum. 
If we choose m^ — 0, we have 



V ^ lml4>' + ^(^^ + n^r + J^i^' + -')0, 




Jf" ^ (j)df'7: - Trdf'(l) + \/ — f 9^7r. (A30) 



Now the potential explicitly breaks the symmetry, and the current is different from the canonical one. The charge Q 
does not commute with the Hamiltonian, and 

This corresponds to full canonical theory after a shift of a to remove its vacuum expectation value. Thus we see that 
the power counting rules allow us to reconstruct the theory without reference to the zero modes. 

Thus in the theory with zero modes dropped and the zero pion mass, the Hamiltonian explicitly breaks the symmetry, 
the vacuum is trivial, and there is no longer the notion of spontaneous symmetry breaking. Current conservation is 
preserved and vastly reduces the number of free parameters present in the effective Hamiltonian constructed from 
power counting and forces the pion to remain massless. 

Of course, we should emphasize that the sigma model we have discussed is only at the tree level. For a complete 
understanding, we need to consider radiative corrections and the subsequent renormalization effects. 

APPENDIX B: MORE ABOUT CHIRAL SYMMETRY ON THE LIGHT FRONT 

In this appendix, we show that the normal chiral transformation on the quark field is inconsistent with the light-front 
constraint equation, where the latter is a result of the light-front equations of motion. We also show that chirality is 
the same as helicity on the light-front. 

Recall that in light-front coordinates, the quark field can be separated into the plus- and minus-components: 
tp ~ tp^ + ?/;-, where tp± — 21^1 i'^ ^^^d only ^_|_ is a dynamical variable. The minus-component V'- is determined 
by the light-front constraint 

V-- = -^(a-L-(«a_L+5^±)+7°mF)V'-H- (Bl) 

The normal chiral transformation is defined by 

ip — > ■!/; + 5^: with 5^ = —iOj^ip. (B2) 

However, in light-front coordinates, only ip^ is a dynamical variable, so that the chiral transformation acts only on 

■0+ — > V'+ + Stp+ , <5')/'+ = -i0j5tp+, (B3) 

and the transformation on ip^ is given by the equation of constraint 

= -iO-fh-77-ai^ ■ (id± + gA±)-tJj+ - i9mFj^j5-KT'^+- (B4) 

10+ ' 10+ 
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Thus, the chiral transformation on tj) on the hght-front is 



= -i6-i5^+ - i9j5T—^a± ■ {id± + gA±)-tJj+ - i9mFj^j5-^4'+- 



id+ 



If we naively use the form S^p — —iOj^ip — — 16*75 (-0+ + V'-), 






9A±)ip+ ~ i6'r7ii?757°---V'+. 



This is obviously inconsistent with (B5), which is a consequence of the light-front equation of constraint. 
For the massive quark theory, 

d^,J5 = 2imFV'75V', 
where ji^ — '07'^ 75?/! is the axial vector current. The axial vector charge on the light-front is then 



Qlf ^ dx tfx_LJ^{x). 



(B5) 
(B6) 

(B7) 
(B8) 



Explicit calculation using the field expansions and normal ordering leads to 



fLF 



dk+cPki 
2(27r)3 



^A[6l(fc)&^(fc)+4(fc)d,(fc)]. 



(B9) 



Thus Q\p measures the helicity. 

In the canonical theory the chiral current is conserved for zero fermion mass. In the theory with the renormalized 
effective Hamiltonian we expect to observe the consequences of spontaneous symmetry breaking, signals of which 
are a conserved axial current and a zero mass pion. However, we do not expect to generate a massless pion unless 
the coupling g exactly matches its renormalized value gs ■ Hence the issues of a massless pion and a conserved axial 
vector current (which we expect to be non-canonical in form) can be addressed meaningfully only after we reach the 
relativistic limit of the theory. 



APPENDIX C: NONRELATIVISTIC LIMIT OF MASSLESS GLUON EXCHANGE 



In this Appendix we discuss the nonrelativistic limit of one gluon exchange; however, we insist on maintaining the 
separate dimensional assignments for longitudinal and transverse coordinates when taking this limit, and we arrive at 
results that differ substantially from other analyses which allow longitudinal and transverse dimensions to mix. Our 
entire program is based on separate power-counting analyses for these dimensions, and if they mix our power counting 
fails. 

In second-order perturbation theory, the effective interaction Hamiltonian in the qq sector due to massless gluon 
exchange is 



where 



H 



(2) 



H 



(2) 
Icrij 1 



H 



(2) 
I(7ij2 



r(2) 






Hpt-ptr 



H 



(2) 
Ia-ij2~ 



-9 ^Q3Qi^Q4Q2-^2y-4- 

Pl 



.y Q3Q1 Q, 



Pi 






Q3Ql-^a4a2^lS3"S2S4 2 / + _ j,+ \2 






1 1 



2y + + o 

Pi -Pi ^ 



m^p+Pl 



''"f+P3± 



1 



iPl±-P3±)^ 
Pt-Pt 



mJ,+P2 

~pr 



{Pl±-P3±V 
Pt-pt 



(CI) 



(C2) 



(C3) 
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and 



M 



2ij- X 



1-1 f'l 

,Pi - Ps 



Pi - Ps ^Pi P3 



Pi Pi 



Xsi 



X~S2 



n Pl -Pi 

+ -f 

Pi ~P:^ 



/ 1 
imp 



1 



Pt Pt 



Pi Pi 



X-s 



(C4) 



Let us denote the total momenta by {P'^,P±). For simplicity, we consider the case P^ — 0, that is, p2i_ = 
—pi±, P4± = —P3±- To get the form of the effective Hamiltonian in the nonrelativistic limit, as an example, we 
consider the matrix element for si — S2 ^ s^ = S4 =|. Then 



Mo 



iPi 



P3±) 



n Pl±iPt-pt) q PIApJ 



Pi 



:>+\2 



24C^p{pf 



\Pi -pI) PiP2\Pi -pI) Pi Pi [Pi -Ps) PiPlPsPi 

where k^ — k^ — ik"^ and k^ = fc^ + ik"^ . This gives, after simplification. 



(C5) 






1 / 4m^ 



2 \pip^ 



. , ,ip-r 



l/4rr4 



Pi P2 Pa Pi 

(p+)2 



+.+^M)- 



1 

El 



Pi P2 P3 Pi 



-p"p"h 



E1E2 V 



Pi, 



ptpt 



pI± 

PsPi 



(C6) 



where 



El = {pii_ - P3iy + 



m. 



-^ipt-ptr iPt-Pt)r-^ ^ 
P1P3 ^ Pi P3 



pI± pI± 



(C7) 



and 



E2 = iPl± - P3A 



^iPt-Ptr-ipt-pt)^^-^^ 
P2 Pi ^ Pi P2 



(C8) 



In the nonrelativistic limit, the term involving PiP^ is negligible compared to the mass term; in both Ei and E2 

'rriF \2 ^„^ "i^F _^ {mF\2 



the third term is negligible compared to the first two terms; furthermore, ^iji^ -^ (^^)'^ and 7^+ ~^ {^^)'^. Thus 

" " Pi P2 P4 ^ 



pfpf 



in the nonrelativistic limit, the effective interaction in the qq sector due to massless one-gluon exchange reduces to 



n_,NR_ 

n-2i3 - 



2rparpa ' 



-Ag^T'^T 



/mF\2 




1 




1 C^^S" 


P3±y + 


i^^npt- 
1 


pt)' 


^Pt^ {P2^ 


-Pi±r 


+ i'fnpt 


-pt)\ 



(C9) 



The Coulomb interaction has the form -^ in momentum space. If we further make the approximation for (C9) that 



P]^ w mp + p\ and pj w mp + p| , we arrive at the familiar Coulomb interaction in momentum space 

1 



Ucrj = -Ag^T^T" 



{Pl - P3f 



(CIO) 



with {pi — ps)^ = {pii_ —P3_Lf' + {Pi —pi)''- The Coulomb potential in position space is then obtained by the Fourier 
transformation. 



(P 



1-P3) Att J r 



(Cll) 
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with r — \x — x'\ the separation between the two fermions. Obviously, He = — J (P x(P x' j^ {x) ^^J"" {x') has lost boost 
invariance but has recovered rotational invariance. 

We wish to represent the interaction 7i^^ in a form where longitudinal boost invariance is kept. To do so, we use 
an interpolating Fourier transformation: 

2 1 



/mF\ 

"+ ^ i{{^ipt-pt)){^V')-iPi±-P3±)-V±} 



i j "it^-yy^^"'^' 



X 



/mp\^ 1 



where y± = xn_ — X2i^ = 5x^_ is the transverse separation of the two constituents, but the longitudinal separation y^ = 
5x is more complicated since the above expression is not a complete Fourier transformation to light-front coordinate 

space. But for quahtative purposes we may define a hght-front "radial" coordinate ri — y {6x±)^ + {^^Y{5x~Y, 

and we see from (C12) that iJ;^^ has the form of the light-front Coulomb potential (2.6) which we have introduced 
in the artificial potential in Section II. E. This light-front radial coordinate is invariant under longitudinal boosts. 
However, it is not invariant under transverse boosts. For more physical remarks, see Section II. E. 

Relativistically, the spinor matrix elements are different in different helicity sectors. In the nonrelativistic limit this 
helicity dependence vanishes, and we get the same interaction in all helicity sectors. 

APPENDIX D: PHASE SPACE CELL DIVISION: EQUAL-TIME CASE 

Let us try to motivate why a phase space cell analysis is useful for the study of Hamiltonians. Recall free field 
theory in equal-time coordinates. A free field is diagonalized in terms of momentum eigenstates. Thus it is customary 
to express free fields in terms of creation and annihilation operators for particles in plane wave states. The free 
Hamiltonian is 

Ho= (fkoJk alak, (Dl) 



with uok = y'lJp + k^. 

Next consider interactions. For illustrative purposes, we will choose a 0^-type interaction, 

H,nt^ f d^x f dPki f d^k2 f d^k3a\ki)a{k2)a{k3)e-'^''^+'''-''''>-''. (D2) 



Before trying to solve the Hamiltonian quantitatively, it is important to make qualitative estimates of various terms 
in the Hamiltonian. For this purpose, the creation and annihilation operators aj, and aj, are quite inappropriate. For 
example, the operator aj. creates a particle in a plane wave state and hence creates a state of infinite norm. Thus 
when we try to estimate the "order of magnitude" of a[. we get infinity! What can be done to avoid this catastrophe? 
It has been suggested |38J (also see [KM ) some time ago that the creation and annihilation operators which depend on 
a continuous variable k be expanded in terms of a discrete, complete, orthonormal set of functions Ui{k). Specifically, 



Ofc 



^^Ui{k)ai , (D3) 



4 = E<w«I- (D4) 



The coefficients of this expansion are themselves creation and annihilation operators for particles in normalizable 
"wave packet" states. Now 
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Ha = ^ Ci,i^ al^Ui^, (D5) 

iii2 

where 

a,,, = f (fk LJk < {k)u,, (fc) . (D6) 

Hf) is now the Hamihonian of an infinite number of osciUators coupled to each other (no longer diagonal). Now 
making Hq diagonal means making self-interactions of the individual oscillators more important than interactions 
between different oscillators. Consider the matrix C. If the u;'s are orthogonal, C would be diagonal except for 
the factor ojk- Even with the factor uok, one should be able to keep the off-diagonal elements small. If the u^'s are 
"properly chosen" (that is, as localized in momentum space as possible) then for distinct momentum shells ii and 
12, the functions Ui-^ and Ui^ do not overlap very much. If Ui-^ and Ui^ are in the same momentum shell but different 
spatial cells, this fact will be reflected in a rapid phase variation of u* Ui^ as a function of fc, which again makes the 
integral small. Thus with a properly chosen set of u^'s, for order of magnitude estimates we need to consider only 
diagonal elements Ci^i^ of C. 

The interaction Hamiltonian in terms of wave packet states is 

Hint = ^ aj^ai^aig Alibis (D7) 

with 

Au.«3= / d^^ J d^k, J d^k2 I d^ksul {k,)u,, (k2)u,, (fcg) e-*('=2+'=^-'=i)- 

d^xv*^{x)vi2{x)vi^{x). (D8) 



If the u's are as localized as possible in momentum space, then the w's are as spread out as possible in coordinate space 
and hence the off-diagonal elements of D are as big as the diagonal elements of D. Thus if we want to approximate 
D by its diagonal elements for order of magnitude estimates, we need w's to be as localized in coordinate space as 
possible. Thus we clearly have to satisfy almost mutually exclusive requirements, which means that we have to settle 
for a compromise. 

A qualitative implementation of the compromise is as follows. We choose Ui{k) such that if Aki is the momentum 
width of the function Ui(k) and Axi is the width of the Fourier transform Vi{x) of Ui{k), then the product AkiAxi is 
near the lower limit set by the uncertainty principle, that is, 

Afc^Ax, « - . (D9) 

The set of functions {ui} is chosen to be complete and orthonormal. What does this mean? Think of Ui as occupying 
a cell of unit volume in phase space. Completeness means that the total volume occupied by the u^'s must fill all 
space. Orthogonality means that the regions occupied by different u^'s do not overlap. 

How are we going to divide the phase space into cells? Consider one dimensional example for visualization. Divide 
the momentum space into an infinite number of cells, the «*'* cell being 

2-' <k< 2''+^ . (DIO) 

For each momentum cell separately, divide the position space linearly into cells of the appropriate size. The position 
space coordinate is chosen to be I < 2~*x < I + 1. 

Let us give some examples. Consider z = 0. Then 1 < fc < 2, and the length of the momentum cell is Lk = 1. The 
position cell division is / < a; < / -f 1. Thus the position space coordinates arc, for example, 

l = -l -l<a;<0, (DU) 

1 = 0< a; < 1, 
; = +1 1< a; < 2, 

and so on, so that each position space cell has length L^ = 1- Thus the " volume" of each phase space cell is LkL^ — 1. 
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Consider i — I. Then ^ < fc < 1, and the length of the momentum cell is Lk = 5. The position cell division is 
I < ^x < I + 1. Thus the position space coordinates are, for example, 

l = -l -2<a;<0, (D12) 

1 = 0< a; < 2, 
l = +l 2<a;<4, 

and so on, so that each position space cell has length L^ = 2. Thus the "volume" of each phase space cell is again 
LkLx = I- 

One can label each phase space cell by {I, i). An illustration of the division of phase space into cells is given in Fig. 
13. 

Why have we chosen a logarithmic scale for momentum cells? One would like to have each momentum cell correspond 
to a distinct energy scale so that order of magnitude estimates become meaningful. A linear choice obviously fails 
for this purpose, while a logarithmic choice is the most obvious one that fulfills it. Note that there is nothing sacred 
about the factor of 2; for a specific example where the factor of 2 appears as the optimum choice, see Ref. [p9[. 
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